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Abstract 

A  dynamic  model  which  simulates  snowmelt  in  Colorado  sub- 
alpine  watersheds  for  all  combinations  of  aspect,  slope,  elevation, 
and  forest  cover  composition  and  density  is  described.  The  model 
simulates  winter  snow  accumulation,  the  energy  balance,  snowpack 
condition,  and  resultant  melt  in  time  and  space.  Detailed  flow 
chart  descriptions  of  the  various  components  of  the  model  and  a 
program  listing  are  presented. 
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Computer  Simulation  of 
Snowmelt  Within  a  Colorado 
Subalpine  Watershed 

Charles  F.  Leaf  and  Glen  E.  Brink 


The  systems  approach  is  a  way  to 
determine  the  probable  effects  of  land 
management  on  the  many  interdependent 
hydrologic  components  in  the  Colorado 
subalpine  forest.  Accordingly,  we  are 
developing  mathematical  models  to  simulate 
the  hydrology  of  this  high-elevation 
ecosystem. 

As  a  first  step,  we  are  modeling:  (1) 
winter  snow  accumulation,  (2)  the  energy 
balance,  (3)  snowpack  condition,  and  (4) 
resultant  melt  in  time  and  space  under  a 
variety  of  conditions.  Combinations  of  aspect, 
slope,  elevation,  and  forest  cover  composition 
and  density  are  included.  The  computer 
program  described  was  initially  written  by  the 
Watershed  Systems  Development  Unit  at  the 
Pacific  Southwest  Forest  and  Range 
Experiment  Station  (Willen  et  al.  1971).  We 
have  revised  the  original  program  to  better 
represent  conditions  in  the  Rocky  Mountain 
region.  With  this  snowmelt  model,  we  have 
simulated  the  probable  effects  of  forest  cover 
manipulation  and  additions  to  the  winter 
snowpack  through  weather  modification. 

The  model  consists  of  three  parts:  (1)  the 
determination  of  the  form  of  precipitation  (rain 
or  snow),  (2)  the  melting  process,  and  (3) 
snowpack  condition  in  terms  of  energy  level 
and  free  water  requirements.  Shortwave  and 
longwave  radiation  represent  the  energy 
available  for  snowmelt.  In  the  forested 
environment,  shortwave  radiation  reaching 
the  pack  is  estimated  by  means  of  a  trans- 
missivity  coefficient  function,  which  depends 
on  the  density  and  composition  of  the  forest 
cover  (Miller  1959,  Reifsnyder  and  Lull  1965). 
Radiation  inputs  are  adjusted  for  slope  and 
aspect  (Frank  and  Lee  1966).  Reflectivity  of  the 
snowpack  is  varied  according  to  precipitation, 
the  energy  balance,  and  time  (U.S.  Army  1956). 

The  snowpack  is  assumed  to  behave  as  a 
dynamic  heat  reservoir;  thus  all  elements  in 
the  snowmelt  portion  of  the  model  are 
expressed  in  units  of  heat.  The  net  external 
energy  balance  is  computed  at  the  snow 


surface.  Rain  and  snow  are  converted  from 
inches  at  the  prevailing  air  temperature  to 
equivalent  gram-calories.  Each  precipitation 
event  is  added  algebraically  as  a  caloric-heat 
event  to  develop  the  heat  reservoir  or 
snowpack.  Temperatures  within  the  snowpack 
are  computed  using  unsteady  heat  flow  theory. 
The  snowpack  will  yield  melt  water  only  when 
it  has  reached  a  zero  energy  deficit  (snowpack 
temperature  =  0°  C  )  and  its  free  water  holding 
capacity  is  satisfied.  Snowmelt  rates  after  the 
pack  is  primed  are  governed  primarily  by  the 
longwave  and  shortwave  energy  balances  at 
the  snow  surface.  The  discussion  which 
follows  is  a  detailed  flow  chart  description  of 
the  various  components  of  the  model.  The 
model  has  been  programed  for  the  CDC  6400 
computer-  at  Colorado  State  University. 


Program  Description 

Program  MELTMOD  (Main  Program) 

The  main  program  (fig.  1)  makes  no 
computations,  but  serves  only  to  coordinate 
the  flow  of  the  program.  Within  the  listing  of 
the  main  program  is  included  a  dictionary  of 
all  variables  that  are  available  for  use  and 
cross  referencing  by  the  subroutines  as  well  as 
the  main  program.  (This  allocation  of 
computer  memory  for  mutual  access  by  many 
routines  is  designated  "blank  common.") 
Other  dictionaries  appear  throughout  the 
subroutines  for  variables  used  exclusively 
within  those  routines.  Blank  common  is 
arranged  so  that  additional  variables  may  be 
added  easily;  by  using  a  very  extensive  blank 
common,  the  routines  were  broken  down  into 
simpler  logical  units  rather  than  one  or  two 


-The  use  of  trade  and  company  names  is  for  the 
benefit  of  the  reader;  such  use  does  not  constitute  an 
official  endorsement  or  approval  of  any  service  or 
product  by  the  U.  S.  Department  of  Agriculture  to  the 
exclusion  of  others  that  may  be  suitable. 
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Figure  1.         1-  Subroutine  AFFECTS  (fig.  2) 

The  decisions  necessary  to  compute  the 
response  to  input  data  are  started  within  this 
routine  and  continued  in  other  routines.  The 
primary  decision  to  be  made  here  is  the  classifi- 
cation of  a  precipitation  event.  An  event  is 
classified  as  follows: 

a.  If  the  daily  minimum  temperature  is  less 
than  or  equal  to  32°  F  ,  or  if  the  daily  maxi- 
mum temperature  is  less  than  or  equal  to 
35°  F  ,  it  is  a  snow  event. 

b.  If  the  minimum  temperature  is  greater  than 
35°  F  ,  it  is  a  rain  event. 

__,  c.  If  the  minimum  temperature  is  between 

32°  F  and  35°  F,  and  the  maximum 
temperature  is  greater  than  35°  F  ,  it  is 
a  mixture  of  rain  and  snow  event.  The  total 
precipitation  is  partitioned  into  rain  and 
snow  by  the  formula  presented  in  subroutine 
MIXTURE. 

The  precipitation  is  added  to  the 
snowpack,  and  the  caloric  input  or  loss  due  to 
the  precipitation  is  calculated.  Effects  on  the 
snowpack    are    computed    by  subroutines 


large  routines.  The  model  contains  10 
subroutines,  discussed  in  the  following  order: 


1.  AFFECTS 

2.  CALIN 

3.  CALOSS 

4.  DIFMOD 

5.  GETREF 


6.  LINK 

7.  MIXTURE 

8.  RADBAL 

9.  RAINED 
10.  SNOWED 


The  parameters  which  describe  the 
initial  conditions  in  the  model  are  solar 
radiation  transmissivity  coefficient,  forest 
cover  density,  an  initial  pack  temperature  and 
water  equivalent,  and  a  threshold  value  for  use 
in  the  reflectivity  subroutine,  GETREF.  The 
data  items  supplied  on  a  daily  basis  for  the 
simulation  are  incoming  shortwave  radiation, 
daily  maximum  and  minimum  temperatures, 
observed  water  equivalent  (optional), 
precipitation,  and  density  of  the  snowpack 
(optional).  Observed  snowpack  temperatures 
for  use  by  the  diffusion  model,  subroutine 
DIFMOD,  are  also  optional,  but  recommended. 


Figure  2. 
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CALIN  or  CALOSS.  Subroutines  GETREF 
and  RADBAL  are  then  called  to  compute  the 
reflectivity  of  the  pack,  and  the  influences  of 
the  radiation  balance  as  computed  from  air 
temperature. 

2.  Subroutine  CALIN  (fig.  3) 

This  routine  computes  the  effects  of 
caloric  input  from  either  the  radiation  balance 
or  from  rainfall.  A  check  is  made  to  see  if  the 
input  satisfies  an  exisiting  calorie  deficit, 
given  by  the  equation 


CALDEF  =  203.2  W„ 


[1] 


where 


CALDEF  =  heat  deficit  of  the  snowpack  in 
calories,  and 

=  the  "cold  content"  which  repre- 
sents the  heat  required  to  raise 
the  snowpack  to  0°  C  ,  in 
equivalent  inches  of  water  (U.S. 
Army  1960). 


^    CALIN  ^ 


Figure  3. 
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If  the  input  energy  is  not  sufficient  to 
make  up  the  heat  deficit,  the  only  action  taken 
is  to  adjust  the  pack  temperature  to 
compensate  for  the  input.  If  the  calorie  deficit 
is  satisfied  and  there  is  still  sufficient  input  to 
generate  melt,  the  melt  is  added  to  the  free 
water  content  of  the  snowpack  and  deleted 
from  the  ice  content.  A  new  free  water  holding 
capacity  is  then  computed  (4  percent  by  weight 
of  the  ice  content)  and  any  excess  free  water  is 
deleted  from  the  pack. 


3.  Subroutine  CALOSS  (fig.  4) 

This  routine  computes  the  effects  of  a 
loss  of  energy  from  the  pack.  The  loss  may  be 
due  to  either  a  negative  radiation  balance  (the 
back  radiation  exceeding  incoming  radiation)' 
or  "cold  content"  added  by  precipitation  (fresh 
snow).  The  snowpack  is  first  examined  to  see  if 
any  free  water  is  present  in  the  pack;  if  not,  the 
calorie  deficit  is  increased  to  compensate  for 
the  loss.  If  free  water  is  present,  the  loss  is  used 
to  freeze  all  or  part  of  it  (80  calories  per  cm.  of 
free  water).  If  the  loss  is  more  than  sufficient  to 
freeze  all  of  the  free  water  the  remaining  heat 
loss  is  used  to  create  a  calorie  deficit  (eq.  1). 
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Figure  4. 
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4.  Subroutine  DIFMOD  (fig.  5) 

This  routine  controls  the  snowpack 
temperature  through  the  winter  months. 
Because  the  snowpack  is  typically  below  0°  C 
during  much  of  the  snow  accumulation  season, 
heat  flow  theory  is  used  to  index  its  thermal 
regime. 

Total  heat  energy  input  to  the  snowpack 
is  represented  by  air  temperature  measured  3 
to  4  feet  above  the  snow  surface.  The 
temperature  in  this  boundary  layer  region  is 
assumed  to  be  an  integration  of  all  the 
processes  involved,  including  incoming  and 
outgoing  radiation,  wind,  temperature,  and 
humidity  in  the  overlying  airmass.  The  heat 
flow  equation  as  solved  by  this  routine  is  given 
by  Quick  (1967)  and  Riley  et  al.  (1969): 


3^ 


=  c_„p 


3Ts 
9t 


[2] 


where 


kg  =  thermal  conductivity  in  cal/°C/cm/sec, 
Cg  =  specific  heat  in  cal/gm/°C, 
Pg  =  snowpack  density  in  gm/cm^ 
Tg  =  snowpack  temperature  in  °C  , 

z  =  depth  within  snowpack  (from  surface) 
in  cm.,  and 

t  =  time  in  sec. 


Equation  [2]  can  be  rewritten  as: 


=  K. 


V 


3Ts 
3t 


[3] 


where 

=  thermal  diffusivity  in  cm^/sec. 

The  thermal  diffusivity  is  assumed  to  vary 
with  density  according  to  the  relation  proposed 
by  Schwerdtfeger  (1963): 


2 

=  (3pi  -  ps)c. 


[4] 


where 


kj^  =  thermal  conductivity  of  ice  in 

cal/°C/cm/sec., 
Pi  =  density  of  ice  in  gm/cm^,  and 
c  -j^  =  specific  heat  of  ice  in  cal/gm/°C. 


Figure  5. 
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The  density  is  assumed  to  be  constant 
throughout  the  pack,  although  it  may  vary 
with  time. 

A  finite  difference  solution  for  equation 
[3]  is  obtained  as  follows  (Smith  1965, 
Richtmyer  and  Morton  1967): 

Let 


3u  9^u 
3l=  ^  9^ 


a  =  const.  >0 


[5] 


where  equation  [5]  is  the  nondimensionalform 
of  equation  [3].  (It  should  be  noted  that  the 
number  representing  the  depth  of  the  pack  is 

The  forward  difference  approximation 
to  equation  [5]  is  given  by 

[6] 

u.  -  u.    .       u.,T    ,-2u.    ,  +  u.   -,  . 


m 


where 


X  =  ih,  (i  =  0,  1,  2,  .  .  .) 
t  =  jm  (j  =  0,  1,  2,  .  .  .) 
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Equation  [6]  can  be  written  as: 

[7] 

u.    .,T=u.    .+TiiK   (u,   -,    .-2u.    .+u.,,  .) 

1,:       V     1-1, J       i,J  x+l,J 

which  gives  the  expHcit  solution  for  the 
unknown  temperature  uij+i  at  the  (i,j+l)  th 
mesh  point  in  terms  of  known  temperatures 
along  the  j  th  time  row.  Hence,  it  becomes 
possible  to  calculate  the  unknown  pivotal 
values  of  u  along  the  first  time  row,  t  =  m(fig.  6) 
in  terms  of  known  boundary  and  initial  values 
along  t  =  0,  then  the  unknown  pivotal  values 
along  the  second  time  row  in  terms  of  the 
calculated  pivotal  values  along  the  first,  and  so 
forth. 

Equation  [7]  utilizes  the  average  air 
temperature  (assumed  to  be  the  temperature  of 
the  surface  of  the  snowpack),  ground 
temperature,  and  temperature  at  the  midpoint 
of  the  snowpack  during  the  previous  interval  to 


simulate  present  snowpack  temperature.  The 
midpoint  temperature  is  defined  initially  from 
observed  data  and  simulated  thereafter.  The 
diffusion  model  has  been  found  to  be  mathe- 
matically stable  only  for  reasonably  deep 
packs  (>4.7  inches  of  water  equivalent).  To 
insure  stability  in  the  model,  each  24-hour 
period  is  divided  into  two  12-hour  intervals  and 
the  result  is  averaged.  The  average  air 
temperatures  for  each  12-hour  interval  are  the 
"high  average"  (the  mean  of  the  maximum 
and  daily  mean  temperatures)  and  the  "low 
average"  (the  mean  of  the  minimum  and  daily 
mean  temperatures.) 

Once  the  snowpack  becomes  isothermal, 
subroutine  RADBAL  computes  the  energy 
balance.  Subroutine  LINK  determines  when 
the  diffusion  model  is  to  be  used  and  when  the 
radiation  balance  is  to  be  used.  The  diffusion 
model  may  be  initialized  or  readied  by  external 
controls  or  when  the  radiation  balance  creates 
a  calorie  deficit. 


t 


Figure  6. 
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5.  Subroutine  GETREF  (fig.  7) 

This  routine  selects  the  reflectivity  for 
use  in  the  radiation  balance  routine,  RADBAL. 
In  essence,  there  are  two  functions  from  which 
the  reflectivity  is  selected:  one  for  the 
accumulation  phase  of  the  snowpack  and  one 
for  the  melt  phase  (U.S.  Army  1960).  The 
reflectivity  is  assumed  to  start  at  a  higher 
value  and  to  decrease  more  slowly  with  time 
during  accumulation  than  during  the  melt 
season.  To  determine  which  function  to  use,  a 
check  is  made  on  the  following  conditions: 

a.  If  no  new  snow  was  added  to  the  pack  dur- 
ing a  given  day,  the  length  of  time  since 
fresh  snow  is  increased  by  1  day  and  the 
function  used  on  the  previous  day  is  used 
again.  The  accumulation  function  computes 
reflectivity  for  15  consecutive  days.  After  15 
days,  control  shifts  to  the  next  lower  value 
in  the  melt  function,  which  computes  reflec- 
tivity for  an  additional  26  days.  After  this 
period  of  time,  the  reflectivity  is  a  constant 
40  percent.  The  minimum  value  for  the  re- 
flectivity in  any  case  is  assumed  to  be 
40  percent. 

b.  If  new  snow  has  fallen  on  the  pack,  a  further 
check  is  made  to  see  if  it  was  a  "spring"  or 
"winter"  snow.  This  is  determined  by  com- 
paring the  maximum  daily  temperature 
against  a  temperature  threshold  which 
varies  according  to  elevation,  forest  cover, 
and  aspect.  If  it  is  warmer  than  the  thresh- 
old, no  change  is  made  and  the  control 
flows  as  explained  in  part  a.  If  it  is  colder, 
indicating  "winter"  snow  conditions,  a- 
nother  check  determines  if  the  snowpack 
is  isothermal  and  thus  selects  which  func- 
tion to  use  while  re-initializing  the  day 
counter. 


6.  Subroutine  LINK  (fig.  8) 

This  routine  is  the  interface  between  the 
thermal  diffusivity  model  and  the  radiation 
balance  routine.  As  discussed  above,  the 
diffusion  model  is  used  to  control  the  snowpack 
temperature  during  the  winter  season.  Air 
temperature  is  the  driving  variable  in  both 
routines,  so  they  are  somewhat  related.  The 
decision  as  to  which  method  to  use  is  basically 
dependent  on  whether  the  snowpack  is  still 
under    "winter"    conditions    or  "spring" 


^    GETREF  J 


Figure  7. 
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conditions.  The  first  check  is  to  see  if  the 
radiation  balance  has  simulated  a  caloric  gain 
or  loss  from  the  pack.  A  loss  requires  a  check  to 
see  if  the  pack  has  been  isothermal  with  free 
water  content  ("spring")  or  if  it  has  been  cold 
("winter").  If  it  was  already  cold,  the  diffusion 
model  is  used  without  further  question.  If  the 
pack  did  contain  free  water,  the  loss  is  used  to 
freeze  it  and,  if  a  calorie  deficit  would  be 
created,  the  diffusion  model  is  re-initialized  to 
isothermal  conditions  and  the  radiation 
balance  is  adjusted  accordingly.  In  the  event  of 
a  caloric  gain,  a  check  is  made  to  see  if  the 
calorie  deficit  would  be  satisfied;  if  not,  the 
diffusion  model  is  used.  If  the  deficit  is 
satisfied,  a  further  check  is  made  to  see  if  at 
least  0.05  inch  of  free  water  is  generated;  if  not, 
the  diffusion  model  is  set  to  isothermal 
conditions  and  control  passes  to  the  diffusion 
model.  However,  if  sufficient  free  water  is 
generated,  the  diffusion  model  is  turned  off 
and  the  radiation  balance  resumes  full  control. 
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Figure  8. 
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7.  Subroutine  MIXTURE  (fig.  9) 


Q  MIXTURE  J       Figure  9. 
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This  routine  partitions  rain  and  snow 
precipitation  components  which  fall  on  days 
when  the  maximum  temperature  is  greater 
than  35°  F  and  the  minimum  temperature  is 
greater  than  32°  F  ,  but  less  than  35°  F.  The 
partitioning  is  performed  through  the 
following  equation: 

SNOW  =  (PRECIP)  (1  -  B/A)  [8] 

where 

B  =  the  difference  between  the  maximum 

temperature  and  35°, 
A  =  the  difference  between  the  maximum 

and  minimum  temperatures. 

The  "average"  temperature  for  the  snow  period 
is  computed  as  the  mean  of  the  minimum 
temperature  and  35°,  while  that  for  rain  is  the 
mean  for  the  maximum  and  35°.  Control  then 


8.  Subroutine  RADBAL  (fig.  10) 

This  routine  uses  air  temperature  to 
compute  the  net  radiation  balance  through  the 
Stefan-Boltzmann  function  as  follows: 


where 

Lp  =  potential  longwave  radiation  at  a 
given  temperature, 
O  =  Stefan-Boltzmann  constant  for  a 
24-hour  period: 

1.17  X  10    (langleys)/(day)/(°K)  -\  and 
T  =  temperature  in  °K. 

The  shortwave  radiation  component  in  the 
radiation  balance  is  computed  as  a  simple 
function  of  the  transmissivity  coefficient  and 
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the  cover  density.  Solar  radiation  is  adjusted 
for  slope  and  aspect,  according  to  tables 
published  by  Frank  and  Lee  (1966).  The 
downward  longwave  component  is  computed 
by  using  the  average  air  temperature  in  the 
Stefan-Boltzmann  function  according  to  the 
equations: 


From  sky  to  snow: 


-7, 


[10] 


=  0^(1  -  Cp)(1.17  X  10    )(T  ^) 


where 


Cj-j=  forest   cover  density  expressed  as  a 
decimal, 

=  ambient  air  temperature  in  °K,  and 

a  =  a  factor  (1  or  0.75)  which  accounts  for 
clear  or  cloudy  skies. 


From  forest  cover  to  snow: 

Lp  =  1.17  X  IQ-'c^T/ 


[11] 


where 

Tp=  radiation  temperature  of  foliage  in  °K. 


Figure  10. 
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The  upward  component  (back  radiation)  is 
computed  by  using  either  the  average  daily  air 
temperature  ("winter"  conditions)  or  the 
minimum  air  temperature  ("spring" 
conditions)  for  the  radiation  temperature  of  the 
snowpack.  In  the  case  of  the  back  radiation, 
however,  the  temperature  used  must  be  less 
than  or  equal  to  0°  C.  The  equations  used  are: 


From  snowpack  to  forest: 


=  1.17  X  10-' 


[12] 


where 


To  =  radiation  temperature  of  the  snowpack  in 
°K. 


From  snowpack  to  sky: 


[13] 


-7. 


Lgg  =  (1  -  C^)(1.17  X  10  'XTg'^) 


Once  the  upward  and  downward 
components  have  been  calculated,  they  are 
combined  to  get  a  net  longwave  balance  as 
follows:  if  the  skies  are  clear  (no  precipitation), 
only  75  percent  of  the  downward  longwave 
radiation  given  by  equation  [10]  is  combined 
with  the  back  radiation  from  the  snowpack 
(U.S.  Army  1960).  The  radiation  balance  under 
the  forest  canopy  is  computed  by  equations  [11] 
and  [12].  During  "winter"  conditions,  Tp  =  Tg 
and  Lf=  LsF-  If  there  was  precipitation,  a 
check  is  made  to  see  if  it  was  snow.  When  there 
is  snow,  the  longwave  balance  is  assumed  to  be 
zero.  Otherwise,  under  cloudy  skies,  the 
downward  component  and  the  back  radiation 
are  merely  combined  algebraically. 

Once  the  net  radiation  is  computed,  the 
decision  is  made  by  the  linking  routine  to  use 
either  the  diffusion  model  or  the  radiation 
balance.  The  effects  of  the  decision  are  then 
computed  by  either  Subroutine  CALIN  or 
Subroutine  (I!ALOSS,  depending  upon  the  net 
gain  or  loss  registered  by  the  radiation 
balance. 


9.  Subroutine  RAINED  (fig.  11) 

This  routine  computes  the  effects  of  a 
rain  event  on  the  snowpack  according  to  the 
equation: 


=  (V  ("^>  <^r) 


[14] 
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Cr 


=  calorie  gain  due  to  rainfall 

=  specific  heat  of  water:  1  cal./gm/°  C  , 


A  T  =  difference  between  temperature  (°C  ) 

at  which  rain  falls  and  0°  C  ,  and 
Pj^  =  amount  of  rainfall  in  centimeters. 

If  the  pack  is  cold,  the  caloric  input  from 
the  rain  is  used  to  satisfy  all  or  part  of  the 
calorie  deficit.  If  the  input  more  than  satisfies 
the  deficit,  the  remainder  is  contributed  as  free 
water  and  the  caloric  input  from  that 
remaining  is  allowed  to  generate  other  melt.  If 
the  pack  was  already  isothermal,  the  entire 
amount  of  rain  is  added  to  the  pack  as  free 
water,  and  the  calories  contribute  to  the  melt 
rate. 


Figure  11. 
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10.  Subroutine  SNOWED  (fig.  12) 

This  routine  computes  the  effects  of  a 
snow  event  on  the  snowpack  according  to  the 
equation 


=  (CgXADCPg) 


where 

L  J  =  calorie  gain  or  loss  due  to  snowfall, 

Cg  =  specific  heat  of  ice:  0.5  cal./gm/°  C  , 

AT  =  difference  between  temperature  (°  C  ) 
at  which  snow  falls  and  0°  C  ,  and 

Pg  =  waterequivalentof  snowfall  in  centi- 
meters. 

If  the  snow  falls  within  the  "warm" 
range  of  32°  to  35°  F ,  no  action  is  taken 
concerning  the  caloric  loss.  However,  snow 
falling  at  lower  temperatures  increases  the 
calorie  deficit,  and  this  change  is  accounted 
for. 


(  SNOWED    ^  Figure  12. 
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[15] 


Unlisted  Routines 

The  routines  listed  below  are  not 
included  in  the  flow-charted  descriptions  since 
they  are  merely  utility  routines  for  handling 
input,  output,  or  internal  processing  of  the 
information.  They  are  not  part  of  the  model; 
rather,  they  are  the  necessary  routines  to 
implement  the  model  on  a  digital  computer. 

INITIAL        RDPACK  STORE 

PLOTTER      READER  WRITER 

A  complete  listing  of  the  snowmelt  model  as 
described  above  is  included  in  appendix  1. 
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Model  Verification 

The  model  has  been  tested  on  field  data 
obtained  from  the  667-acre  Deadhorse  Creek 
watershed  (fig.  13)  at  the  Fraser  Experimental 
Forest  (Leaf  1971).  Snowmelt  rates  at  eight 
locations  on  Deadhorse  Creek  were  re- 
constituted (fig.  14).  Agreement  between 
observed  and  simulated  melt  rates  was  good  at 
all  locations.  The  model  is  structured  so  that  a 
minimum  number  of  variables  must  be  ad- 
justed to  obtain  a  satisfactory  fit.  In  this  case, 
only  two  variables  were  adjusted  to  obtain 
satisfactory  agreement:  (1)  cover  density  (a 
vegetation  type  and  composition  parameter), 
and  (2)  a  shortwave  radiation  transmissivity 
coefficient.  Figure  15  compares  observed  and 
simulated  snowpack  accumulation  and  melt  at 
the  high-elevation  north  slope  and  low- 
elevation  south  slope  sites  on  Deadhorse  Creek 
in  1969.  Table  1  is  an  example  of  the  computer 
output. 


With  the  cover  density  and  trans- 
missivity coefficients  fixed  as  determined  in 
1969,  the  model  has  given  good  results  in 
simulating  snowmelt  on  Deadhorse  Creek  for 
the  1964-71  runoff  seasons.  Although  it  is  a 
simplification  of  the  real  life  system,  the  model 
produces  reliable  results. 

Empirical  studies  have  shown  that 
various  watershed  management  practices 
exert  a  significant  effect  on  snowmelt  rates 
and  resultant  streamflow.  Because  the  model 
described  here  is  a  mechanistic  representation 
of  the  snowmelt  process,  its  careful  use  should 
enable  the  resource  manager  to  better  under- 
stand how  a  given  management  alternative 
will  affect  snowmelt  before  it  is  implemented. 
We  believe  it  to  be  a  useful  tool  for  predicting 
the  probable  effects  of  land  management 
practices  on  the  timing  and  amount  of 
snowmelt. 


Figure  13.  —  Aerial  view 
(foreground)  of  Deadhorse 
Creek  watershed,  Fraser 
Experiment  Forest.  The 
north-facing  slope  is  in  the 
lower  lefthand  corner  of 
the  photograph. 
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Figure  14.  —  Locations  where  field  measurements  were  taken  on  Deadhorse  Creek. 
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Figure  15.  —  Observed  and  simulated  snowmelt  on  upper  north  (A)  and  lower  south  (B) 
slopes  during  1969  snowmelt  runoff  season. 
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TABLE  1 


TRANSMISSIVITY  COEFFICIENT    =  .25 


TEMPERATURE 

(Fl 

DATE 

MAX 

MIN 

AVE 

4 

7 

69 

30 

0 

12 

0 

21.0 

4 

8 

69 

32 

0 

13 

0 

22 . 5 

4 

9 

69 

36 

0 

18 

0 

27.  0 

4 

10 

69 

40 

0 

23 

0 

31  .  5 

4 

11 

69 

45 

0 

27 

0 

36.  0 

4 

12 

69 

34 

0 

25 

0 

29.  5 

4 

13 

69 

33 

0 

23 

0 

28.  0 

4 

14 

69 

36 

0 

22 

0 

29.  0 

4 

15 

69 

40 

0 

20 

0 

30.0 

4 

16 

69 

28 

0 

15 

0 

21  .  5 

4 

17 

69 

32 

0 

15 

0 

23.  5 

4 

18 

69 

40 

0 

16 

0 

28.  0 

4 

19 

69 

44 

0 

22 

0 

33 . 0 

4 

20 

69 

47 

0 

17 

0 

32.  0 

4 

21 

69 

53 

0 

30 

0 

41.5 

4 

22 

69 

45 

0 

33 

0 

39.  0 

4 

23 

69 

55 

0 

32 

0 

43.5 

4 

24 

69 

52 

0 

24 

0 

38.0 

4 

25 

69 

24 

0 

9 

0 

16.5 

4 

26 

69 

20 

0 

8 

0 

14.0 

4 

27 

69 

27 

0 

8 

0 

17.5 

4 

28 

69 

38 

0 

15 

0 

26.  5 

4 

29 

69 

42 

0 

22 

0 

32.  0 

4 

30 

69 

47 

0 

26 

0 

36.  5 

ELEVATION  10500  FT. 

ASPECT  NE,  SLOPE  35  PERCENT 


SNOWMELT  RUNOFF  SIMULATION  MODEL 
FRASER  EXPERIMENTAL  FOREST,  COLORADO 
DEADHORSE  CREEK,  UPPER  NORTH  SLOPE 

COVER  DENSITY 

PRECIP    (IN)         NET  RAD    (CAE)         ENERGY  SNOWPACK  PREDICTED 


DAY 

ACCOM 

SHORT 

LONO 

BAL 

(CAE) 

TEMP  (c) 

W.E 

(I 

0.11 

13 

45 

8.3 

-7.3 

0 

1 

-1.0 

13. 

45 

0.00 

13 

45 

13.  5 

-18.0 

-4 

5 

-1.3 

13 

45 

0.00 

13 

45 

29.4 

-49.  5 

-20 

1 

-2.4 

13 

45 

0.  00 

13 

45 

30.7 

-21.3 

9 

5 

-1.9 

13 

45 

0.  00 

13 

45 

25.7 

-20.6 

5 

1 

-1.6 

13 

45 

.29 

13 

74 

8.9 

-5.6 

2 

8 

-1.4 

13 

74 

0.  00 

13 

74 

28.4 

-22.4 

6 

1 

-1.0 

13 

74 

0.00 

13 

74 

31  .4 

-30.8 

6 

-1.0 

13 

74 

.  76 

14 

50 

24.2 

-30.8 

-7 

7 

-1.4 

14 

50 

.11 

14 

61 

3.1 

-.7 

1 

5 

-1.3 

14 

61 

0.  00 

14 

61 

22.9 

-22.4 

6 

-1.3 

14 

61 

0.  00 

14 

61 

23.7 

-24.7 

-1 

0 

-1.3 

14 

61 

0.  00 

14 

61 

34.3 

-32.0 

2 

3 

-1.2 

14 

61 

0.  00 

14 

61 

38.6 

-36.6 

2 

0 

-1.1 

14 

61 

0.  00 

14 

61 

40.  3 

-38.0 

2 

4 

-1.0 

14 

61 

0.00 

14 

61 

34.1 

-22.5 

11 

6 

-.3 

14 

61 

.52 

1  5 

13 

37.4 

-31  . 1 

6 

2 

0.  0 

15 

13 

0.00 

15 

13 

45.7 

-1.0 

44 

7 

0.0 

15 

13 

0.00 

15 

13 

36.5 

-27.4 

9 

1 

0.0 

15 

13 

.14 

15 

27 

16.9 

0.0 

15 

1 

0.0 

15 

27 

0.00 

15 

27 

33.0 

-18.5 

14 

5 

0.0 

15 

27 

0.00 

15 

27 

45.0 

-11.2 

33 

7 

0.0 

15 

27 

0.00 

15 

27 

48.3 

-19.8 

28 

5 

0.0 

15 

13 

0.00 

15 

27 

56.6 

-18.  5 

38 

1 

0.0 

14 

94 
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Appendix  I:  Complete  Listing  of  Snowmelt  Model 

Program  MELTMOD  c  between  the  surface  and  the  ground,  and  location  3  IS 

c  the  ground  temperature. 

kROGRAM  MELTMDD   I  input, output, TAPE5  =  INPUT,TAPE5=0UTPUT I  C            SIMTEM2   =  THE   SNOWPACK  TEMPERATURE   AT  THE  MIDDLE  NODE   FOR  THIS 

C  THIS   IS   A  REWRITTEN  VERSION  OF   -SNMELT-,    A  SNOW   ACCUMULATION  MOL/EL  C  INTERVAL,    AS  SIMULATED  BY   SUBROUTINE  DIFMOD 

C               DEVELOPED  BY  THE   U.    S.   FOREST   SERVICE   AT   BERKELEY.      THE   RADIATION  C  SIHTEH3  =   THE   AVERAGE  SNOWPACK  TEMPERATURE   FOR   THIS    INTERVAL,  AS 

C               BALANCE,   AS  DERIVED  FROM  AIR  TEMPERATURE,    SERVES   AS   THE  MODEL  C  DERIVED  BY  SUBROUTINE  DIFMOD 

C               PARAMETERS   AND  CHARACTERISTICS.  C  SNOMELT   -  MELT  OELIVEO   IN   INCHES  FOR  THE  INTERVAL 

C   C  SOBSEOV  -  ARRAY  FOR  STORING  THE   OBSERVED  WATER   EQUIVALENT  FOP 

C  DICTIONARY  OF   BLANK  COMMON  C  PLOTTING 

C   C  SPRECIP  -  ARRAY  FOR   STORING  THE   PRECIP  FOR  PLOTTING 

C  ACTDATE   -  THE   DATE   OF  THE   RECORDING  OF   THE   VALUES   IN  -ACTUAL-  C            SPREOV  -  ARRAY  FOR   STORING   THE   PREDICTED  WATER   EQUIVALENT  FOR 

C  ACTUAL   -   THE   OBSERVED   SNOWPACK   TEMPERATURES  USED  TO   INITIALIZE   THE  C            SUBTITL  -  ONE   CARD  SUBTITLE,   SIMILAR  TO  -TITLE- 

C                                       DIFFUSION  MODEL    (SUBROUTINE  DIFMDDI.      LOCATIONS   1  C  TCQEFF  -  THE   TRANSMI SI  V I T Y  COEFFICIENT   USED  TO  ESTIMATE  THE  NET 

C                                       THROUGH   19  ARE   THE   OBSERVED  TEMPERATURES   AND  LOCATION  C  SHORT  WAVE   RADIATION  REACHING  THE   SNOWPACK.  SEE 

C                                       20   IS   THE   GROUND  TEMPERATURE    (-XX. X   DEGREES   C).  C  REIFSNYDER  AND  LULL,   RADIANT  ENERGY   IN  RELATION  TO 

C                                       LOCATION  21   IS  THE   DISTANCE   BETWEEN   THE   MEASURED  C  FORESTS,   USES   TECH.   BUL   13**,  1965. 

C                                       SNOWPACK  TEMPERATURES   (XX. X   INCHES   -  USUALLY   6    IM.)  C  TEMPMAX  -  THE   MAXIMUM  TEMPERATURE  DURING  THE   INTERVAL    IN  DEGREES 

C  AVETEMC  -   THE   DEGREES  CENTIGRADE   EQUIVALENT   OF  -AVETEMF-  C  FARENHEIT 

C  AVETEMF  -  MEAN  OR   AVERAGE  OF  THE  MAXIMUM  AND  MINIMUM  TEMPERATURE  C            TEMPHIN  -   THE   MINIMUM  TEMPERATURE  DURING  THE   INTERVAL    IN  DEGREES 

C                                       IN  DEGREES  FARENHEIT  C  FARENHEIT 

C  BASTEMF   -   BASE   TEMPERATURE  DEGREES   FARENHEIT,   RAIM   TURNS   TO   SNOW  C           THRSHLD  -  THE   THRESHOLD  TEMPERATURE   FOR  DETERMINING  WHETHER  OR  NOT 

C  CALAIR  -  POTENTIAL   LONGWAVE  CALORIC    INPUT   AT   AIR   TEMPERATURE  C                                       TO  RE-INITIALIZE  THE  REFLECTIVITY   FUNCTION  WHEN 

C  CALDEF  -  THE  CALORIE  DEFICIT   IS   THE   NUMBER   OF  CALORIES  NEEDED  C                                       THERE   IS  A  SNOW   EVENT.      IF  THE   MAXIMUM  TEMPERATURE  IS 

C                                       TO  BRING  THE   SNOWPACK  TEMPERATURE  TO  ZERO  DEGREES  C  GREATER  THAN   THE   THRESHOLD  VALUE  DO  NOT  RE'I N I T I AL I ZE 

C                                       CENTIGRADE    (NOTE  SHOULD  BE  MADE   THAT   IT   IS  A  POSITIVE  C  THE   FUNCTION  REGARDLESS  OF  THE  PRECIPITATION 

C                                       QUANTITY)  C  TITLE  -  ONE  CARD  TITLE    (IF  THE   INFORMATION   IS  CENTERED  ON  THE  CARD 

C  CALORIE  -  CALORIES   OF  HEAT   ABSORBED  OR  RELEASED  BY   THE   SNOWPACK  c                                       IT   WILL  BE   PROPERLY  CENTERED  ON  THE  PAGE) 

C                                       FROM  THE  NET   RADIATION   BALANCE  C  TOTPREC   -   THE   ACCUMULATED  TOTAL   PRECIPITATION   IN  INCHES 

C  CALSNOW  -   POTENTIAL   LONGWAVE  CALORIC   LOSS   AT  SNOW   TEMPERATURE  C            USEMEAN   =  0,   USE   MAXIMUM  AND  MINIMUM  TEMPERATURES  AS  READ 

C  CDVDEN  -  THE   COVER   DENSITY   IS  THE  FRACTION  OF  THE  GROUND  OR   SNOW  C                            =1,   REPLACE   THE   MAXIMUM  AND  MINIMUM  TEMPERATURES  BY  THEIR 

C                                       SURFACE   SHADED  FROM  DIRECT   SUNLIGHT  OR  RADIATION  C  MEAN 

C  DATE  -   THE   DATE   BEING  PROCESSED   IN  MMDDYY   FORMAT  C            XMAX  -   MAXIMUM  OBSERVED  OR  PREDICTED  WATER  EQUIVALENT,   USED  FOR 

C  DATES  -  AN  ARRAY  OF   THE   DATE   BEING  PROCESSED.      THE   FIRST  WORD   IS  C  SCALING 

C  THE   MONTH,   THE   SECOND  THE   DAY,    AND  THE   THIRD   THE  YEAR  C  

C  DEN  -  THE  SNOWPACK   DENSITY  READ  FROM   INPUT   CARDS  COMMON  ACTDATE , ACTUAL ( 2 1 ), AVETEMC, AVETEMF 

C  DENSITY  -   THE   DENSITY  OF   THE   SNOWPACK  USED   IN  THE   DIFFUSION  MODEL.  COMMON  BASTEMF 

C                                       IF  -DEN-   IS  ZERO  OR   BLANK,   -DENSITY-   IS  COMPUTED  AS  A  COMMON  C ALA  I R , CALDEF ,CALORI E , CALSNOW , COVDEN 

C                                       FUNCTION  OF   THE   PREDICTED  WATER  EQUIVALENT  COMMON  DATE , DATE S ( 3 ) , DEN, DENSI TY, DREADY 

C  DREADY  =  0,   DIFFUSION  MODEL    (SUBROUTINE   DIFMOD)   NOT    INITIALIZED  COMMON  ENGBAL 

C                          =1,   DIFFUSION  MODEL   INITIALIZED  AND  READY  FOR   SNOWPACK  COMMON  FOQTNOT ( 16 ) , FREEWAT 

C                                       TEMPERATURE   SIMULATION  COMMON  HOLDCAP 

C  ENGBAL  -  THE   TOTAL  CALORIC   INPUT  TO  OR  LOSS   FROM  THE   SNOWPACK  COMMON   I  DATE ( 372 ), I  SNOW , I  TABLE 

C                                       DURING  AN   INTERVAL.      IT   IS  THE   ALGEBRAIC   SUM  OF   THE  COMMON  KOUNT 

C                                       ENERGY   INVOLVED  WITH  THE   PRECIPITATION  AND  THAT   OF  COMMON  LASTUSD , L I NES 

C                                       THE   RADIATION  BALANCE,   -  CALORIE-  COMMON  NEXTACT 

C  FOOTNOT  -  ARRAY  OF  FOOTNOTES  TO  BE   PRINTED  AT  THE   BOTTOM  OF   EACH  COMMON  OBSWEQV 

C                                       PAGE.      TWO  CARDS   ARE  READ,   THE   FIRST   130  CHARACTERS  COMMON  PACKTEM , PARTICE , PAST  I  NT , PLOTOBS , PLOTWE , PREC I P, PREWEQV 

C                                       FORMING  ONE  LINE   AND  THE  LAST   30  CHARACTERS  FORMING  A  COMMON   RAD i N , RADLWN, RAOSWN , REFLECT 

C                                       SECOND  LINE  COMMON   S I MTEHl ( 3 ) , S I MTEH2 , S I HTEM3, SNOMELT , SOBSEQV ( 372 ) , 

C  FREEWAT  -  THE   FREE   WATER  BEING  HELD   BY  THE   SNOWPACK  1                SPREC I P ( 372 ) , SPREQV ( 372 ) , SUBTI TL ( 8 ) 

C  HOLDCAP  -  THE   FREE  WATER  HOLDING  CAPACITY   OF   THE   SNOWPACK  COMMON   TCOEFF , TEMPMAX , TEHPHI N, THRSHLD, T ITLE ( 8 ), TOTPREC 

C                                       (ASSUMED  TO  BE   FOUR  PERCENT  OF   THE   WATER   EQUIVLAENT)  COMMON  USEMEAN 

C  IDATE  -  ARRAY  FOR   STORING  THE   DATES   FOR   PLOTTING  COMMON  XMAX 

C  ISNOW  -  A   SWITCH  TURNED  ON  BY  SUBROUTINE   SNOWED  WHEN  THE   PRECIP  INTEGER  ACTDATE 

C                                       WAS  SNOW  AND  THEN  OFF   BY   SUBROUTINE  GETREF  AFTER  INTEGER   DATE , DATES , DREADY 

C                                       COMPUTING   THE   REFLECTIVITY  FOR  THE  GIVEN   INTERVAL  INTEGER  FOOTNOT 

C  ITABLE   =   0,   NO  PRINTING  OF  TABULATED  RESULTS  FROM  THE   SIMULATION  INTEGER   PASTI NT, PLOTOBS , PLOTWE 

C                          =1,   PRINT  THE   TABULATED  RESULTS  FROM  THE   SIMULATION  INTEGER  SUBTITL 

C  KOUNT  -  COUNTER  FOR  THE   NUMBER   OF  CARDS   READ,   MAXIMUM  OF   372   DUE  INTEGER  TITLE 

C                                       TO  THE   DIMENSIONS   OF  THE  VARIABLES  FOR   STORING   THE  INTEGER  USEMEAN 

C                                       INFORMATION  FOR   PLOTTING  COHMON/CONVERT/F I VE9TH, THI RTY2 

C  LASTUSD  -  AN   INDICATOR   USED   IN  SUBROUTINE   GETREF  TO  DETERMINE  DATA  F I VE9TH, THIRTY2/. 5555555556 , 32 . 0/ 

C  WHICH  REFLECTIVITY   FUNCTION  TO  USE  C  INITIALIZE   THE   MODEL  AND  READ  THE  PARAMETER  CARDS 

C  LINES  -   THE  LINE  COUNTER  FOR  PAGE   EJECTION  10  CALL  INITIAL 

C  NEXTACT  -  A  ONE  COLUMN  CODE   (COL  70)   ON  THE   INPUT  CARDS   TO  C  READ  A   DATA  CARD 

C                                       INDICATE  WHEN  A  CARD  CONTAINING  THE   ACTUAL   SNOWPACK  20  CALL  READER  (lEND) 

C  TEMPERATURES   IS   TO   BE   READ.      WHEN  COLUMN  70   IS  NOT  C  A  BLANK  CARD  MAY  BE  USED  TO  SEPARATE   SETS  OF  DATA 

C                                       BLANK  OR   ZERO,   CONTROL   SHIFTS   TO  SUBROUTINE  RDPACK  I F ( I  END. NE . 0. OR .DATE .LE . 0  I    GO  TO  60 

C  WHICH  THEN  READS  THE  NEXT  CARD  C  SEE  HOW   THIS   INTERVAL  AFFECTS   THE  SIMULATION 

C  OBSWEQV  -  OBSERVED  WATER  EQUIVALENT   OF   THE   SNOWPACK   IN   INCHES  CALL  AFFECTS 

C  PACKTEM  -   THE   EFFECTIVE   TEMPERATURE   OF  THE   SNOWPACK  C  IF   THE   TABLE   IS  BEING  PRINTED,   WRITE  THIS  LINE 

C  PARTICE  -  THE   PORTION  OF   THE   PREDICTED  WATER  EQUIVALENT   THAT   IS  IF(ITABLE)  *0,'»0,30 

C                                       ICE.      THIS  QUANTITY  PLUS  FREE   WATER   IS   THE   TOTAL  30  CALL  WRITER 

C  PREDICTED  WATER  EQUIVALENT    (-PREWEQV-)  C  IF   THE   PLOT   IS  TO  BE   DONE,    STORE  THIS    INFORMATION.      THEN  GO  ON  TO 

C  PASTINT  -  NUMBER  OF    INTERVALS  SINCE   THE  LAST   INITIALIZATION  OF   THE       C   THE   NEXT  CARD 

C                                       REFLECTIVITY  FUNCTION  'lO   IF(PLOTWE)  20,20,50 

C  PLOTOBS  =  0,   00  NOT   PLOT  THE   OBSERVED  WATER   EQUIVALENT  50  CALL  STORE 

C                                 1,    PLOT  THE  OBSERVED  WATER  EQUIVALENT   (OPERATIVE   ONLY   IF  GO  TO  20 

C  -PLOTWE-   IS  TURNED  ON)  C  ALL  CARDS  HAVE   BEEN  READ,    SO  PLOT  THE  SIMULATION  AFTER  WRITING 

C  PLOTWE   =  0,   DO  NOT   PLOT  THE   SIMULATION  C   THE   FOOTNOTES  ON  THE   LAST  PAGE 

C                          =1,   PLOT  THE   SIMULATION,    PRECIP,   ETC.  60   I F ( I  TABLE. NE. 0 )    WRITE    (6,910)  FOOTNOT 

C  PRECIP  -  OBSERVED   PRECIPITATION   IN   INCHES  -ilO  FORMAT  (  IH013A10/ 1X3A 10 ) 

I  PREWEQV  -  PREDICTED  WATER   EQUIVALENT   OF  THE   SNOWPACK    IN   INCHES  I F ( PLOTWE . NE . 0 )   CALL  PLOTTER 

C  RAOIN  -  RADIATION   IN   IS   THE   TOTAL   INCIDENT   SHORT  WAVE  RADIATION  C  IF   THE   END   OF   FILE  HAS   NOT   BEEN  SENSED,   GO  ON  TO  THt  NEXT  SET  OF 

C  RADLWN  -  NET  LONG   WAVE   RADIATION   IS   THE   ALGEBRAIC   SUM  OF  THE  LONG         C   DATA 

C                                       WAVE   RADIATION  FROM  THE   FOREST   AND  THE  LONG  WAVE  IF(IEND)  70,10,70 

C                                       RADIATION  LOST  BY  THE   SNOWPACK   TO  THE   CANOPY  70  STOP 

;  RADSWN  -   THE  CALORIC   INPUT   TO  THE   PACK   BY   THE  NET   SHORT  WAVE  END 
C  RADIATION 

C  REFLECT  -   THE   FRACTION  OF   RADIATION   THAT   IS   REFLECTED  BY   THE   SNOW  .  ^ 

C                                       AS  DERIVED  BY  SUBROUTINE   GETREF  bUDrOUtine  ArrlliL/lb 

C  SIMTEMl   -  AN   ARRAY  USED  PRIMARILY   IN   SUBROUTINE   DIFMOD   IN  THE 

C                                       SIMULATION  OF   THE   AVERAGE   SNOWPACK   TEMPERATURE.  SUBROUTINE  Arl-tCTS 

C  TO  INSURE   STABILITY  OF   THE   DIFFUSION  MODEL,    THE  C  DETERMINE   THE   EFFECTS  OF  THE   DATA  FROM  THli  CARD 

C                                       DAY   IS   PARTITIONED   INTO   12  HOUR    INTERVALS,   AS  COMMON   ACTDATE , ACTUAL ( 2 1 ), AVETEMC , AVETEMF 

C                                       DISCUSSED   IN   SUBROUTINE   DIFMOD.      THIS   ARRAY   STORES  COMMON  BASTEMF 

C                                       THE   CONDITIONS   PRESENT   DURING  THIS    INTERVAL   FOR   USE  COMMON  C ALAI R , CALDEF .CALOR I E , CALSNOW, COVDEN 

C                                       IN  THE   SIMULATION  ON  THE  NEXT    INTERVAL.      LOCATION   1  COMMON  DATE , DATE S ( 3 ) ,DEN, DENSI TY , DREADY 

C                                       STORES  THE   AVERAGE   AIR  TEMPERATURE    (ASSUMED  TO  BE  COMMON  ENGBAL 

C                                       THE   SURFACE   TEMPERATURE  OF  Tl^r  i  SNOWPACK  )  ,   LOCATION   2  COMMON  FOOTNOT  (  16 ),  FREEWAT 

C                                       IS  THE   SNOWPACK  TEMPERATURE   *   |  .  NODE   MIDWAY  COMMON  HOLDCAP 
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COMMON  IDATE(372) ,ISNOW,ITABLE 

COMMON  KOUNT 

COMMON  LASTUSDiLINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON  PACKTEM,PARTICEiPASTINT,PLOTOBS,PLOTWE,PRECIP,PREHE0V 
COMMON  RAO  I N.RADLHN.RADSWN, REFLECT 

COMMON  SIMTEMlOl  ,  S I MTEM2  ,  S I MTEM3,  SNOMELT,  SOBSEQVI 372  )  , 
1  SPRECIP(372) ,SPRE0V(372) .SUBTITL(8) 

COMMON  TC0EFF,TEMPMAX,TEMPMIN,THRSHLD,TITLE(8),T0TPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER   DATE, DATES, DREADY 
INTEGER  FOQTNOT 

INTEGER  PASTINT,PLOTOBS,PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

COMMON/CONVERT /F I VE9TH,TH I RTY2 

C  START  THE   ENERGY   BALANCE  AND  SNOHMELT   AT  ZERO  FOR  THIS  INTERVAL, 

C   BUT  ACCUMULATE   THE  PRECIPITATION 

ISNOW  =  0 

ENGBAL   =  0.0 

SNOMELT  =  0.0 

TOTPREC   =  TOTPREC  +  PRECIP 

C  IF  A  SNOWPACK  ALREADY  EXISTS,   GO  FIND  THE   EFFECT   OF  THIS  I>4TERVAL 

C  ON  THE  SNOWPACK 

IF(PREWEOV)  I0,10,*0 

C  SINCE   THERE   IS   NO  SNOWPACK,   CHECK  TO  SEE   IF  THIS   INTERVAL   HAS  ANY 

C   PRECIPITATION  AND   IF  SO,    SEE   IF   IT   IS   ALL  SNOW 

10   IF(PRECIP)  120,120,20 

C  IF   THE   MINIMUM  TEMPERATURE    IS   BELOW   FREEZING,    CONSIDER  THE 

C   PRECIPITATION  TO  BE   ALL  SNOW 

20   IF(TEMPMIN.LE.THIRTY2.0R.TEMPMAX.LE.BASTEMF)   GO  TO  60 

C  SEE   IF   THE   PRECIPITATION   IS   ALL  RAIN  OR   A  MIXTURE  ON  BARE  GROUND 

IF(TEMPMIN  -  BASTEMF)  90,30,30 

C  IT   IS  ALL  RAIN,   DO  NOT   START   BUILDING  UP  THE  SNOWPACK 

30  SNOMELT  =  PRECIP 
GO  TO  120 

C  A   SNOWPACK  EXISTS.      IF   THERE   IS   PRECIPITATION,    DETERMINE  THE  TYPE. 

C   BUT  OTHERWISE,   JUST  GO  ON  TO  COMPUTE   THE  REFLECTIVITY  AND 

C   RADIATION  BALANCE 

*0   IFIPRECIP)  110,110,50 

C  THERE   IS   PRECIPITATION  ON  AN  EXISTING  PACK.      IF   IT    IS  NOT  ALL 

C   SNOW,   GO  SEE   IF   ANY  OF   IT  WAS  SNOW 

50   IF(TEMPMIN.GT.THIRTY2. AND. TEMPMAX.GT. BASTEMF)   GO  TO  70 

60  CALL   SNOWED   (AMINl    ( AVETEMC , 0.0 ) , PREC I P ) 
GO  TO  110 

C  SEE   WHETHER  THE  PRECIPITATION  ON  AN   EXISTING  PACK  WAS   ALL   RAIN  OR 

C   A  MIXTURE   OF   RAIN  AND  SNOW 

70  IF(TEMPMIN  -  BASTEMF)  90,80.80 

C  THIS   IS  A  RAIN  ON  SNOW   EVENT.      THE  TEMPERATURE   FOR   COMPUTING  THE 

C   DEPLETION  OF   THE  TOTAL  CALORIE   DEFICIT   IS  THE  DIFFERENCE  OF  THE 

C   AVERAGE   TEMPERATURE  AND  FREEZING   (0.0  DEGREES  CENTIGRADE) 

80  CALL  RAINED    ( AVETEMC , PREC I P ) 
GO  TO  100 

C  THIS   IS  A  MIXTURE   OF   RAIN  AND  SNOW  EVENT 

90  CALL  MIXTURE 

C  IF  THE   PACK  WAS   ENTIRELY  MELTED,   BYPASS  COMPUTATION  QF  THE 

C   REFLECTIVITY  AND  THE  RADIATION  BALANCE 

100   IF(PREWEOV)  120,120,110 
C  GET  THE   REFLECTIVITY  FOR  THIS  INTERVAL 

110  CALL  GETREF 

C  COMPUTE   THE  RADIATION  BALANCE  AND   ITS   EFFECT  ON  THE  PACK 

CALL  RADBAL 
RETURN 

C  THERE   IS  NO  SNOWPACK  -  REDEFINE   THE  RADIATION  BALANCE  TO  A 

C   NEGATIVE  VALUE  TO  ASSURE  THE   PROPER  SELECTION  OF  THE  REFLECTIVITY 

C   FUNCTION   IN  SUBROUTINE  GETREF  WHEN  THERE   IS  A  SNOWPACK 

120  CALORIE  =  -1.0 

RE  rURN 

END 


Subroutine  CALIN 


SUBROUTINE   CALIN  (CALORIN) 

C  THIS   SUBROUTINE  COMPUTES  THE  EFFECTS   OF  THE  CALORIC   INPUT  ON  THE 

C   SNOWPACK 

COMMON  ACTDATE, ACTUAL ( 21) , AVETEMC , AVETEMF 

COMMON  BASTEMF 

COMMON  CALAIR,CALDEF,CALDRIE,CALSNDW,COVDEN 
COMMON  DATE,DATES(3) , DEN, DENSI TY , DREADY 
COMMON  ENGBAL 

COMMON  FO0TNOT(16) .FREEWAT 
COMMON  HOLDCAP 

COMMON   IDATE(372) , ISNOW, ITABLE 

COMMON  KOUNT 

COMMON  LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEOV 

COMMON   PACKTEM,PAR  "^ICE  ,  PAST  I  NT ,  PLOTOBS  ,  PLOTWE  ,  PREC  I  P,  PREWEOV 
COMMON  RAD  I N,RADLWN,RADSWN, REFLECT 

COMMON   SIMTEH1I3) , SI MTEM2 , S I MTEM3 , SNOMELT , SOBSEOV ( 372 ) , 
1  SPR£CIP(372) ,SPREQV(372) ,SUBTITL(8) 

COMMON   TCOEFF,TEMPHAX,TEMPMIN,THRSHLD, TITLE  I  8), TOTPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE, DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER  PASTINT,PLOTOBS,PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 


INTEGER  USEMEAN 
C  ADD  THESE   CALORIES   INTO  THE   ENERGY  BALANCE 

ENGBAL   =  ENGBAL   »  CALORIN 
C  SEE   IF   A  CALORIE  DEFICIT  EXISTS    IN  THE  PACK 

COMPARE   =  CALORIN  -  CALDEF 

IF(COMPARE)  10,20,30 

C  THERE   IS   A  CALORIE  DEFICIT,    BUT  THE   INPUT   DID  NOT  COMPLETELY 

C   WIPE   IT  OUT.      ALL   OTHER  CONDITIONS  ARE  UNCHANGED 

10  CALDEF  =  -  COMPARE 
C  1.27  =  0.05  *  2.5A 

PACKTEM  =  C0MPARE/(PREWE0V*1.27) 

RETURN 

C  THE  CALORIE  DEFICIT  WAS  WIPED  OUT,    BUT  ALL  OTHER  CONDITIONS  ARE 

C   UNCHANGED 

20  CALDEF   =  0.0 

PACKTEM  =  0.0 

RETURN 

C  ANY  DEFICIT  WHICH   DID  EXIST  WAS   WIPED  OUT.      COMPUTE  THE  POTENTIAL 

C   MELT  FROM  THE   REMAINING  CALORIES    ( CALOR I ES/ ( 80 .0  *  2.5*)) 

30  POTMELT  =  COMPARE/203.2 

CALDEF   =  0.0 

PACKTEM   =  0.0 

C  IF  THE   INPUT  WAS   ENOUGH  TO  MELT  THE   WHOLE   PACK,   CONTRIBUTE  THE 

C   WATER   EQUIVALENT   TO  THE  SNOWHELT  AND  ZERO   ALL  CONDITIONS 

IF  (  POTMELT. LT.PARTICE)    GO  TO  <,0 

SNOMELT   =   SNOMELT  +  PREWEQV 

PREWEOV  =  0.0 

PARTICE   =  0.0 

FREEWAT  =  0.0 

HOLDCAP   =  0.0 

RETURN 

C  DEPLETE   THE   ICE   PACK  BY  THE   AMOUNT  MELTED  AND  CONTRIBUTE  THAT 

C  AMOUNT  TO   THE   FREE  WATER 

40  PARTICE   =   PARTICE  -  POTMELT 

FREEWAT  =  FREEWAT  +  POTMELT 

C  COMPUTE   THE   NEW  HOLDING  CAPACITY  OF  THE  PACK  AND  COMPARE   IT  WITH 

C   THE   FREE   WATER  TO   SEE    IF   SNOWHELT   IS  PRODUCED 

HOLDCAP   =  0.04  *  PARTICE 

COMPARE  =  FREEWAT   -  HOLDCAP 

IF(COMPARE.LE.O.O)  RETURN 

C  THE   SNOWMELT  CONTRIBUTED   IS    IN  -COMPARE-.      REDUCE  THE  FREE  WATER 

C   TO  LEAVE   A   PSIHED   PACK  AND  REDUCE  THE  PREDICTED  WATER  EOUIVLAENT 

PREWEQV  =  PREWEQV  -  COMPARE 

SNOMELT  =  SNOMELT  +  COMPARE 

FREEWAT  =  HOLDCAP 

RETURN 

END 

Subroutine  CALOSS 

SUBROUTINE   CALOSS  (CALOUT) 

C  THIS   SUBROUTINE  COMPUTES  THE   EFFECTS  OF  THE  CALORIC  LOSS  ON  THE 

C   SNOWPACK 

COMMON  ACTDATE, ACTUAL(21) , AVETEMC, AVETEMF 

COMMON  BASTEMF 

COMMON  CALAIR, CALDEF, CALORIE, CALSNOW,COVDEN 
COMMON  DATE,DATES(3) , DEN, DENS ITY , OR EADY 
COMMON  ENGBAL 

COMMON  FOOTNOTdb)  .FREEWAT 
COMMON  HOLDCAP 

COMMON   IDATE(372) . ISNOW. ITABLE 

COMMON  KOUNT 

COMMON   LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON   P AC KTEM.PARTICE.PASTI NT. PLOTOBS. PLOT WE, PRECIP, PREWEQV 
COMMON  R AD  I N.RADLWN.RAOSWN. REFLECT 

COMMON   SIMTEMl (3) . S I MTEH2 . S I MTEM3 , SNOMELT , SOBSEOV I  372 ) , 
1  SPRECIP1372) ,SPREQV(372),SU8TITL18) 

COMMON   TCOEFF.TEMPMAX.TEMPMI N,THRSHLD.TITLE(8) .TOTPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE, DATES , DREADY 
INTEGER  FOOTNOT 

INTEGER   PASTINT, PLOTOBS, PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

C  ADD  ALGEBRAICALLY  THESE  CALORIES   INTO  THE  ENERGY  BALANCE 

ENGBAL  =   ENGBAL  +  CALOUT 

C  SEE   IF   THERE   IS  ANY  FREE  WATER   IN  THE   PACK.      IF  NOT,   THE  LOSS  IS 

C   JUST  CONTRIBUTED   TO  THE  CALORIC  DEFICIT  OF  THE  SNOWPACK. 

C-   REMEMBER   THAT  -CALOUT-   IS  NEGATIVE 

IF(FREEWAT.GT.O.O)   GO  TO  10 

CALDEF  =  CALDEF  -  CALOUT 

GO  TO  50 

C  COMPUTE   THE  CALORIC   LOSS  NECESSARY  TO  FREEZE   ALL   OF   THE   FREE  WATER 

C   (FREE   WATER  *  80.0  *  2.54) 

10  CALNEED  =  FREEWAT  *  203.2 

C  NOW  COMPARE   THAT  NECESSARY   LOSS   WITH  THE  ACTUAL   LOSS.      IF   THEY  ARE 

C   THE  SAME,   THE  FREE  WATER   IS  WIPED  OUT  BUT   NO  OTHER  CONDITIONS  ARE 

C   ALTERED 

COMPARE   =  CALOUT  +  CALNEED 
IF(COMPARE)  20.30.40 

C  THE  LOSS  WAS  MORE   THAN   ENOUGH  TO  FREEZE   IT.     THE -BALANCE  CREATES 

C   AN  ENERGY  DEFICIT   IN  THE   PACK  AND  THE  FREE  HATER   IS  WIPED  OUT 

20  CALDEF  =  -  COMPARE 
30  PARTICE  =  PARTICE   +  FREEWAT 
FREEWAT   =  0.0 
GO  TO  50 

C  ONLY  PART  OF   THE   FREE  WATER  FROZE.     COMPUTE   THE   BALANCE  REMAINING 

C   BALANCE  =   EXISTING  FREE  WATER  -  AMOUNT   FROZEN.  WHERE 

C   AMOUNT  FROZEN  =  CALOR I ES/ ( 80 .0   »  2.54) 

40  FROZEN   =  -  CALDUT/203.2 
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PARTICE  =  PARTICE  +  FROZEN 
FREEWAT  =  FREEWAT  -  FROZEN 
RETURN 

 COMPUTE   THE   NEW  PACK  TEMPERATURE   AND  HOLD  CAPACITY 

50   PACKTEM  =  -CALDEF/ ( PREWEOV*! . 27 ) 
HOLDCAP   =  0.0',  *  PARTICE 
RETURN 
END 


Subroutine  DIFMOD 


SUBROUTINE  DIFMOD 

C  THIS   SUBROUTINE  WAS   DERIVED  FROM  PROGRAM  SIMTEM,   A  SNOWPACK 

C   TEMPERATURE  DIFFUSION  MODEL  DEVELOPED  BY  LEAF    (1970   STUDY  PLAN 

C   FS-RM-1602,   NO.   224,   RMF+RES).      USING  THE   AVERAGE  SURFACE  TEMP 

C   AND  THE   GROUND  TEMP   AS   BOUNDARY  CONDITIONS,   THE  NEW  AVERAGE 

C   SNOWPACK  TEMPERATURE   IS  CALCULATED 

COMMON  ACTDATE , ACTUAL (211 .AVETEHC.AVETEMF 

COMMON  BASTEMF 

COMMON  CALAIR, CALDEF, CALOR I E.CALSNOW.COVOEN 
COMMON  DATE, DATESO)  , DEN, DENSITY, DREADY 
COMMON  ENGBAL 

COMMON  FOOTNCTI 16) .FREEWAT 
COMMON  HOLDCAP 

COMMON  IDATE(372),ISNOW,ITABLE 

COMMON  KOUNT 

COMMON  LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON  PACKTEM, PARTICE, PAST  I  NT, PLOT DBS, PL0TWE,PREC1P,PREWEQV 
COMMON  RA0IN,RADLWN,RA0SWN, REFLECT 

COMMON   SIMTEHl (3) ,SIMTEM2.SIMTEM3,SN0MELT,S0BSEQV( 372), 
1  SPRECIP(372) ,SPRE0V(372) ,SUBTITL(B) 

COMMON  TCOEFF,TEMPMAX,TEMPMIN,THRSHLD,TITLE( 8) .TOTPREC 
COMMON  USEHEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE, DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER  PASTINT,PLOTOBS,PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

C  

C  DICTIONARY 

C 

C  CALDM  -  THE  CALORIC   INPUT  OR  LOSS  AS  DERIVED  BY  THE  DIFFUSION 

C  MODEL 

C  CONSTl   -  THE  FIRST  CONSTANT   IN  THE  EQUATION   FOR   THE  SIMULATION 

C  C0NST2  -   THE   SECOND  CONSTANT   IN  THE  EQUATION  FOR   THE  SIMULATION 

C       _H  -  THE  DISTANCE   BETWEEN  NODES    (CORRESPONDS   TO  THE  -H-    IN  THE 
C  STUDY  PLAN) 

C  

C0HH0N/C0NVERT/FIVE9TH,THIRTY2 

C  IF  THIS    IS   THE   DATE   OF   THE   ACTUAL  OR  OBSERVED  PACK  TEMPERATURES, 

C   BYPASS   THE   SIMULATION  AND  USE   THE  OBSERVED  AVERAGE 

IFIDATE  -  ACTDATE)  10,100,10 

C  IF  THE   DENSITY  WAS  NOT   SPECIFIED  ON   THE   INPUT  CARD,   COMPUTE    IT  AS 

C   A  FUNCTION  OF   THE   HATER  EQUIVALENT.      (THE   FUNCTION  WAS  DERIVED 

C   FROM  OBSERVED  CONDITIONS  ON  THE  ERASER  EXPERIMENTAL  FOREST) 

10  IF(DEN)  20,20,30 

20  DENSITY  =    (EXP((0.0179   *  PREWEOV)   +  3.02))/100.0 
GO  TO  'lO 

30  DENSITY  =  DEN 

C  COMPUTE   THE  DISTANCE   BETWEEN  THE  TWO  NODES    IN  CENTIMETERS 

C  DEPTH  =  PREHEQV/OENSITY 

C  H  =  (DEPTH/2)*2.5'i 

'iO  H  =    (PREWEQV/DENSITY)    *  1.27 

C  THE   THERMAL   DIFFUSIVITY   IS   CALCULATED  FROM  THE  FUNCTION 

C   KV  =   0.01/((2.751  -  DENSITY)*  0.*8).      MATHEMATICAL  STABILITY 

C   REQUIRES   THAT   THE   VALUE  OF   THE  QUANTITY   (INTERVAL    IN   SECONDS  » 

C  KV/H*»2)    BE   LESS   THAN  0.5.      WHEN  A   2't  HOUR    INTERVAL   IS  USED,  THE 

C   SNOW  DEPTH  MUST  EXCEED  30   INCHES   (20   PERCENT  DENSITY)   TO  ACHIEVE 

C   STABILITY.      IN  ORDER  TO   INSURE   STABILITY  WITH  SOMEWHAT  SHALLOWER 

C   PACKS    (ABOUT   18   INCHES),   THE  DAY   IS  DIVIDED   INTO  2  TIME  INTERVALS 

C  OF   12  HOURS   ('V3200  SECONDS) 

C  CONSTl   =   (43200  *  0.0l/((2.751   -  DENSITY)    *  0.48))/H*»2 

CONSTl  =  900. 0/( (2.751   -  DENS  I TY ) *H*H ) 

C  THE   MINIMUM  WATER  EQUIVALENT  WHICH  WILL   ACHIEVE  STABILITY  USING 

C   THE   ABOVE   DENSITY  FUNCTION   IS   4.7  INCHES 

IFICONSTl  -  0.5)  60,50,50 

C  THE   MODEL   IS  UNSTABLE  -   INDICATE   THAT   IT   IS  NOT   READY   FOR  USE  NOW. 

C   (IT  MAY  BE   INITIALIZED  AGAIN  BY   AN  OBSERVED  PACK   TEMPERATURE  CARD 

C   AND   STABILITY  WILL   BE   ASCERTAINED  FROM  THE   WATER   EQUIVALENT  AT 

C  THAT  TIME) 

50  DREADY  =  0 

RETURN 

C  GET   THE   SECOND  CONSTANT 

60  C0NST2   =   1.0  -  CONSTl  -  CONSTl 

C  PERFORM   THE   SIMULATIONI    IN   TWO  PARTS    (ONE  FDR  EACH   12  HOUR  PERIOD). 

C   -SIMTEMl-  HOLDS   THE   THREE   TEMPERATURES   FROM  THE  PREVIOUS  INTERVAL 

C   THAT  ARE  NEEDED  TO  SIMULATE   SIHTEM2,    THE   NODE   AT   THE  CENTER  OF 

C   THE   PACK.      SIMULATE   THE   FIRST   12  HOURS  NOW 

S1MTEM2   =   (CONSTl   ♦    (SIMTEMKl)    +  SIHTEM1(3)))    t   (C0NST2   *  SIMTEMl 
1(2)  ) 

C  THE   AVERAGE   SNOWPACK   TEMPERATURE   IS  THE   AVERAGE  OF   THE   2  NODES 

C   (MIDDLE   AND  GROUND)    IN  BOTH   INTERVALS.      GROUND  TEMPERATURE  IS 

C   CONSTANT,    SO  START   THE   AVERAGE  NOW 

SIMTEM3   =   SIMTEM1(3)   +   SIMTEMl(3)   ♦  SIMTEM2 

C  RESET  -SIMTEMl-   TO  THE   TEMPERATURES   OF  THE    INTERVAL   JUST  SIMULATED 

C   FOR  USE    IN  THE   SECOND   12  HOUR   INTERVAL   SIMULATION.      THE  SURFACE 

C   AIR   TEMPERATURE    IS   SPLIT   INTO  A  LOW  AVERAGE    (    (MEAN+HIN)/2    )  AND 

C   A  HIGH  AVERAGE    (    (HEAN+MAX)/2    )    FOR  USE  WITH  THE  TWELVE  HOUR 

C   INTERVALS.      USE   THE   LOW  AVERAGE  NOW 


SIMTEMKl)    =  AMINl    (  0.  0  ,  (  (  (  TEMPMI  N-TH I  RTY2  )  *F  I VE9TH  )  +  AVET  EMC  )  /  2  .  0  ) 

SIMTEM1(2)    =  SIMTEM2 

C  SIMULATE   THE   SECOND   12   HOURS   AND  COMPUTE  THE   AVERAGE  SNOWPACK 

C   TEMPERATURE 

SIMTEH2   =    (CONSTl   «    (SIMTEMKl)    ♦  SIMTEMK3)))    +   (C0NST2  *  SIMTEMl 
1(2)  ) 

SIMTEM3   =    (SIMTEM3  ♦  SIMTEM2)/4.0 

C  RESET  -SIMTEMl-  USING  THE   HIGH  AVERAGE   FOR  USE  ON  THE  FIRST 

C   INTERVAL   OF  THE   NEXT  DAY 

SIMTEMKl)    =  AMINl    (  0.  0  ,  (  (  (  TEMPMAX-TH I RTY2  )  *F  IVE9TH  )+AVET  EMC  )  /  2  .0  ) 

SIHTEM1(2)    =  S1HTEM2 
C  CHECK   TO  SEE   IF   THE  GROUND  TEMPERATURE   SHOULD  BE  RAISED 

IF(SIHTEM3  +   1.5)  100,80,70 
70   IF(SIMTEM3  +  0.5)  80,90,90 

80   IF(SIHTEMK3).LT.-0.5)    SIMTEM1(3)   =  -0.5 

RETURN 
90  SIMTEM1(3)    =  0.0 
100  RETURN 

END 


Subroutine  GETREF 


SUBROUTINE  GETREF 
C  GET   THE  REFLECTIVITY 

C  

C  DICTIONARY 

C 

C  REFACUM  -  A  REFLECTIVITY  FUNCTION  FOR  THE   SNOWPACK   DURING  THE 

C  ACCUMULATION  PHASE  OF   THE  SNOWPACK 

C  REFMELT  -  A  REFLECTIVITY  FUNCTION  FOR  THE   SNOWPACK   DURING  THE 

C  MELT   PHASE  OF  THE  SNOWPACK 

C  

COMMON  ACTDATE, ACTUAL(21  )  .AVETEMCAVETEMF 
COMMON  BASTEMF 

COMMON  CALAIR, CALDEF, CALOR I E,CALSNOW,COVDEN 
COMMON  DATE,DATES(3) , DEN, DENS  I TY , DREADY 
COMMON  ENGBAL 

COMMON  FOOTNOTI 16) , FREEWAT 
COMMON  HOLDCAP 

COMMON  I0ATE(372),1SN0W,ITABLE 

COMMON  KOUNT 

COMMON   LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON  PACKTEM,PARTICE,PASTINT,PLOTOBS,PLOTWE,PRECIP,PREWEQV 
COMMON  RAD  I N,RADLWN,RADSHN, REFLECT 

COMMON   SIHTEMK3)  ,  S I  MTEM2  ,  S  I MTEM3 ,  SNOMELT  ,  SOBSEQV  (  372  )  , 
1  SPRECIP(372) ,SPREQV(372) ,SUBTITL(B) 

COMMON   TC0EFF,TEMPMAX,TEMPMIN,THRSHLD,TITLE(8) ,TOTPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE, DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER  PASTINT,PLOTOBS,PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

DIMENSION  REFACUM ( 1 5 ), REFMELT ( 15 ) 

DATA  REFACUM/. 80,    .77,    .75,    .72,    .70,    .69,    .68,    .67,    .66,  .65, 
1    .64,    .63,    .62,    .61,  .60/ 

DATA  REFMELT/. 72,    .65,    .60,    .58,    .56,    .54,    .52,    .50,    .48,  .46, 
1   .44,    .43,    .42,    .41,  ;40/ 

C  INCREASE   THE   INTERVAL  COUNTER  BY   I   AND  SEE   IF   THERE   WAS  ANY  SNOW 

PASTINT  =  PASTINT  +  1 
IF(ISNOW)  10,10,80 

C  USE   THE   SAME  FUNCTION  AS  LAST  TIME 

10   IF(LASTUSD)  20,20,50 

C  ACCUMULATION  PHASE   -  AFTER   15   DAYS,   USE  THE   MELT  FUNCTION 

C   STARTING  AT  THE  FOURTH  DAY 

20  IF(PASTINT  -  15)  30,30,40 
30  REFLECT  =   REF ACUM ( PASTINT ) 

RETURN 
40  PASTINT   =  4 
LASTUSD  =  1 
GO  TO  70 

C  MELT  FUNCTION  -  AFTER   15  DAYS,   USE   A  CONSTANT  40  PERCENT 

50  IF(PASTINT  -   15)  70,70,60 
60  PASTINT   =  15 
70  REFLECT   =  REFMELT ( PAST  I  NT ) 
RETURN 

C  THERE    IS  NEW  SNOW  -  DETERMINE   IF   THE   FUNCTION   IS  TO  BE  RE- 

C   INITIALIZED 

80   IF(TEMPMAX  -  THRSHLD)  90,90,10 

C  IT   IS.    SO   SEE   WHICH  FUNCTION   IS  TO  BE  USED 

90  PASTINT   =  0 

IF(PACKTEM)  100,110,110 
100  REFLECT   =  0.91 
LASTUSD   =  0 
RETURN 

C  THE   PACK   IS   ISOTHERMAL,   BUT   IF   THE   ENERGY  BALANCE   FROM  THE 

C   PREVIOUS   INTERVAL   WAS   NEGATIVE,   USE  THE   ACCUMULATION  PHASE 

C   FUNCTION  ANYWAY 

110   IF(CALORIE)  100,120,120 
120  REFLECT   =  0.81 
LASTUSD  =  1 
RETURN 
END 


Subroutine  INITIAL 


SUBROUTINE  INITIAL 
C  READ  THE   PARAMETERS   AND   INITIALIZE  THE  MELT  MODEL 
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COMMON  ACTDATE,ACTUAL(21I , AVETEMC, AVETEMF 
COMMON  BASTEXF 

COMMON  C4LAIR,CALDEF,CAL0RIE,CALSN0W,C0VDEN 
COMMON  DATE .DATES  13) .OEN, DENSITY, DREADY 
COMMON  ENGBAL 

COMMON   F00TN0T(16) iFREEWAT 
COMMON  HOLOCAP 

COMMON   IDATE(372) ,ISNOH,ITABLE 

COMMON  KOUNT 

COMMON   LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  QBSWEOV 

COMMON  PACKTEM,PARTICE.PASTINT,PLOTOBS,PLOTWEtPRECIP.PREHEQV 
COMMaN  RADIN,RADLKN,RADSWN,REFLECT 

COMMON  SIMTEMl (3) . S I MTEM2 . S I MTEM3 , SNOMELT , SOBSEOV I  372 ) . 
1  SPRECIP(372) ,SPRE0V(372) ,SUBTITL( 8) 

COMMON  TC0EFF,TEMPMAX,TEMPMIN,THRSHLD,TITLE(8) .TOTPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE, DATES, DREADY 
INTEGER  FDOTNOT 

INTEGER  PASTINT,PLOTGBS,PLOTHE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

C  ESTABLISH   THE   STANDARD  BASE  TEMPERATURE 

BASTEMF   =  35.0 

C  INITIALIZE   THOSE   VARIABLES  WHICH  ARE  NOT   READ  IN  OR  OTHERWISE 

C   INITIALIZED  BEFORE   BEING  USED 

FREEWAT  =  0.0 

HOLDCAP   =  0.0 

SUM  =  0.0 

TOTAL  =  0.0 

XMAX  =  0.0 

DREADY  =  0 

I  SNOW  =  0 

KOUNT  =  0 

LASTUSD  =  0 

PASTINT  =  0 

LINES  =  999 

C  START  THE   RADIATION  BALANCE  WITH  A  NEGATIVE  VALUE  FOR  POSSIBLE  USE 

C   BY  SUBROUTINE  GETREF   IN  DETERMINING  WHICH  REFLECTIVITY  FUNCTION 

C  TO  USE 

CALORIE  =  -1.0 

C  READ  THE   INSTRUCTION  PARAMETERS   AND  THE  CONTROLS   ON  THE  MODEL 

READ    (5.910)  ITABLE.PLOTWE.PLOTOBS.USEMEAN.TCDEFF.COVDEN.PACKTEM, 
1  PREWEOV.THRSHLD 
910  FDRMAT(«.I1.1X,5F5.2) 
IF(E0F(5))  10,10.20 

C  READ  THE   TITLE.   SUBTITLE  AND  FOOTNOTE  CARDS 

10  READ   15.920)    T I TLE . SUBTITL , FOOTNOT 
920  FORMATIBAIO) 

C  INITIALIZE   THE   ICE   CONTENT   AND  ACCUMULATED  PRECIPITATION 

PARTICE  =  PREWEQV 
TOTPREC   =  PREWEOV 

C  CALCULATE  THE  CALORIE  DEFICIT   FROM  THE  PACK  TEMPERATURE 

C  CALDEF   =  -    (PACKTEM  *  80)   *   (PREWEOV  *  2.54)    /  160 

CALDEF   =  -   PACKTEM  *  PREWEQV   *  1.27 
RETURN 
20  STOP 
END 


Subroutine  LiNK 

SUBROUTINE   LINK  (IRETURN) 

C  THIS   SUBROUTINE   IS  THE   INTERFACE  BtTWEEN  THE  RADIATION  BALANCE 

C   (SUBROUTINE  RADBAL)    AND  THE  DIFFUSION  MODEL    (SUBROUTINE  OIFMOD) 

COMMON  ACTDATE, ACTUAL(21) , AVETEMC, AVETEMF 

COMMON  BASTEMF 

COMMON  CALAIR,CALDEF,CALORIE,CALSNOW,COVDEN 
COMMON  DATE,DATES(3) , OEN, DENS ITY , DREADY 
COMMON  ENGBAL 

COMMON  FOOTNOT ( 16 ), FREEWAT 
COMMON  HOLDCAP 

COMMON   IDATE(372) ,ISNOW,ITABLE 

COMMON  KOUNT 

COMMON  LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON  PACKTEM, PART  ICE, PASTI NT, PLaTOBS,PL0TWE,PRECIP, PREWEOV 
COMMON  RAOIN.RADLWN.RADSWN, REFLECT 

COMMON   SIMTEMl (3) , SI MTEM2 , S I MTEM3, SNOMELT, SOBSEOV( 372 ) , 
1  SPRECIP(372),SPREQV(372) , SUBTITL! 8) 

COMMON  TCOEFF,TEMPMAX,TEMPMIN,THRSHLD,TITLE( 8 ) .TOTPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER   DATE, DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER   PASTINT, PLOTOBS.PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

COMMON/CON VERT/F I VE9TH, THIRTY 2 
C  SEE   IF  THE   RADIATION  BALANCE   IS   AN  ENERGY   LOSS   OR  GAIN 

IF(CALORIE)  10,10,80 

C  THERE   hAS   A  LOSS.      IF  THIS   IS   STILL   WINTER    (NO  FREEWATER),  JUST 

C   GO  AHEAD  AND   USE   THE  DIFFUSION  MODEL 

10   IF(FREEWAT)  20,20,50 

C  USE  THE   DIFFUSION  MODEL  TO  SIMULATE  THE  CURRENT  AVERAGE  SNOWPACK 

C  TEMPERATURE 

20  CALL  OIFMOD 

IF(DREADY)   AO, 40, 30 


C  NOW  MAKE  ANY  NECESSARY  ADJUSTMENTS   IN  THE  RADIATION   BALANCE  TO 

C   CAUSE   THE   PACK  TEMPERATURE  TO   BE  THE  SAME   AS  -SIMTEM3-.     START  BY 

C   GETTING  THE  DIFFERENCE   BETWEEN  THE   LAST   PACK  TEMPERATURE   AND  THIS 

C   ONE,    THEN  CONVERT   IT  TO  CALORIES 

30  CALDM  =    (SIMTEM3  -   PACKTEM)    »  PREWE(3V   *  1.27 

C  ADJUST   THE   LONG  HAVE   PORTION  OF  THE  RADIATION  BALANCE  BY  THE 

C   DIFFERENCE   BETWEEN  THE   CALORIES   DERIVED  FROM  THE  DIFFUSION  MODEL' 

C   AND  THE   ENERGY  BALANCE 

RADLWN   =  RADLWN  +   (CALDM  -  CALORIE) 

CALORIE   =  CALOM 
40   IRETURN  =  0 

RETURN 

C  THE   LOSS   IS   USED  TO  FREEZE   PART   OR  ALL  OF  THE  FREE  WATER,    BUT  IT 

C   MAY  NOT  CREATE  COLD  CONTENT.      IF   IT  WOULD  CREATE  COLD  CONTENT, 

C   RE-INITIALIZE  THE   DIFFUSION  MODEL  TO  0  AND  ADJUST  THE  ENERGY 

C   BALANCE  ACCORDINGLY 

50  CALL  CALOSS  (CALORIE) 

IF(FREEWAT  -  0.5)  50,60,70 
60  SIMTEMKl)    =  AMINl    1  AVETEMC, 0.0 ) 

SIMTEM1(2)    =  0.0 

SIMTEMl(3)    =  0.0 

DREADY  =  1 

C^  MAKE   ANY  NECESSARY  ADJUSTMENTS  TO  THE  ENERGY  BALANCE  TO  COMPENSATE 

C   FDR   THE  COLD  CONTENT  THAT  WOULD  HAVE   BEEN  GENERATED  BY   THIS  LOSS 

C   AND   ZERO  THE   COLO  CONTENT 

CALORIE   =  CALORIE  +  CALDEF 

RADLWN   =  RADLWN  +  CALDEF 

PARTICE   =   PARTICE  +  FREEWAT 

PREWEQV  =   PREWEQV  +  FREEWAT 

FREEWAT   =  O.C 

CALDEF  =  0.0 

PACKTEM  =  0.0 
70   IRETURN  =  1 

RETURN 

C  THERE   IS  CALORIC   INPUT  TO  THE  PACK.     CHECK  TO  SEE   IF  CONDITIONS 

C   INDICATE   THAT  THE  DIFFUSION  MODEL  SHOULD  BE  TURNED  OFF  AND  THE 

C   ENERGY  BALANCE  USED  FOR  SPRINGTIME  SIMULATION.      CONSIDER  FIRST 

C   ANY  COLD  CONTENT   (INCLUDING  THAT  OF  THE   PREVIOUS  DAY  AND  ANY 

C   CREATED  BY  A   SNOW  EVENT  ON  THIS  DAY).      IF  THERE   IS   COLD  CONTENT, 

C   CHECK   THE   AVERAGE  AIR  TEMPERATURE   AND  THE  SNOHPACK  TEMPERATURE 

C   FROM  THE  PREVOUS   DAY  FOR   ARBITRARILY  CHOSEN  SPRINGTIME  CONDITIONS 

C   AND   IF   ALL  ARE   NOT   SATISFIED,   GO  AHEAD  AND  USE  THE  DIFFUSION 

C  MODEL 

80   IF(CALDEF)  170.170.90 

90   IF (AVETEMC. LE. 0.0. OR. PACKTEM. LE. -0.7)    GO  TO  20 

C  SINCE   SPRINGTIME  CONITIDNS   PREVAIL,    RECOMPUTE  THE   BACK  RADIATION 

C   AND   THE  NET  RADIATION   BALANCE    (REMEMBER,    IF  THERE   IS  SNOW,  THE 

C   LONGWAVE   IS  ASSUMED  TO  BE  ZERO,   SO  THERE  WOULD  BE  NO  NEED  TO  MAKE 

C   ANY  ADJUSTMENTS) 

IF(ISNOW)  100,100,140 
100  USE  =    (TEMPMIN  -  THIRTY2]    *  FIVE9TH 
IF(USE.GT.O.C)    USE   =  0.0 

CALSNOW  =   1.17E-7  *    ((USE  +   273.15)    *•  4) 
IF(PRECIP)  110.110.120 
110   RADLWN  =   {(1.0  -  CDVDEN)   »   ((0.757  •  CALAIR)   -  CALSNOW))   +  (COVDEN 
1  •   (CALAIR  -  CALSNOW)) 
GO  TO  130 
120  RADLWN  =  CALAIR  -  CALSNOW 
130  CALORIE   =  RADSWN  *■  RADLWN 

C  RE-INITIALIZE  THE   DIFFUSION  MODEL  TO  THESE  CONDITIONS   (BUT   IF  THE 

C   INPUT   IS  MORE   THAN  ENOUGH  TO  WIPE  OUT  THE  CALORIE   DEFICIT.  JUST 

C   LET   IT  BRING   THE   PACK  TQ   ISOTHERMAL.      IN  THIS  HAY,    TWO  CONSECU- 

C   TIVE  DAYS  OF   INPUT  ARE  REQUIRED  TO  GENERATE  FREE  WATER) 

140  ACTDATE   =  DATE 

COMPARE  =  CALORIE  -  CALDEF 
IF(COMPARE)  160,150.150 

C  INITIALIZE   THE   DIFFUSION  MODEL  TO   ISOTHERMAL  CONDITIONS 

150   SIMTEMKl)   =  0.0 
SIMTEM1(2)   =  0.0 
SIMTEM1(3)   =  0.0 
SIMTEM3  =  0.0 
GO  TO  20 

C  REDEFINE  THE  SURFACE  TEMPERATURE  AND  COMPUTE  THE  NEW  AVERAGE  PACK 

C   TEMPERATURE.      THEN  COMPUTE  THE   MIDDLE  NODE   AS   A  FUNCTION  OF  THAT 

C   AVERAGE,   THE   SURFACE  TEMPERATURE   AND  THE  GROUND  TEMPERATURE 

C   (WHICH  REMAINED  UNCHANGED) 

150  SIMTEMKl)    =  AHINl    ( 0.  0  ,  AVETEMC  ) 

SIMTEM3  =  COMPARE/ ( PREWEOV   *  1.27) 

SIMTEM1(2)    =    (3.0  »  SIMTEH3)   -  SIMTEMKl)   -  SIHTEM1(3) 
SIMTEM1(3)   =  0.0 

GO  TO  20 

C  THERE   IS  INPUT  TQ  THE  PACK  AND  THE  PACK   IS   ALREADY   ISOTHERMAL.  IF 

C   THIS   ENERGY  WILL  CREATE  AT  LEAST  0.05   INCH   (ARBITRARY  AMOUNT)  OF 

C   FREE  HATER,   TURN   THE  DIFFUSION  MODEL  OFF   AND  LET  THE  ENERGY 

C  BALANCE   TAKE   ITS  COURSE 

170   IFIFREEHAT   +   (CALORIE/203.2)   -  0.05)  150,180.180 
180  DREADY  =  0 
IRETURN   =  0 
RETURN 
END 


Subroutine  MIXTURE 


SUBROUTINE  MIXTURE 
— THIS   SUBROUTINE  CONTROLS  THE  COMPUTATIONS  FOR   A  PRECIPITATION 
  EVENT   THAT   IS  A  MIXTURE  OF  SNOW  AND  RAIN 

 DICTIONARY 

AMTSNOW  -  THE  AMOUNT  OF   PRECIPITATION  OCCURRING  AS   SNOW  (INCHES) 
TFORAIN  -  THE  TEMPERATURE  FOR  COMPUTING  THE   DEPLETION  OF  THE  TOTAL 

CALORIE  DEFICIT  CAUSED  BY  THE  RAIN   (DEGREES  C) 
TFORSNO  -   THE   TEMPERATURE   FOR  COMPUTING  THE  CONTRIBUTION  OF  THE 

SNOW  TO  THE  TOTAL  CALORIE  DEFICIT    (DEGREES  C) 
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c  

COMMON  ACTDATE .ACTUAL (21) , AVETEMC , AVETEMF 
COMMON  BASTEMF 

COMMON  CALAIR,CALDEF,CALORIE,CALSNOW,CQVOEN 
COMMON  DATE.DATESO)  ,  DEN,  DENS  I TY  ,  DR  EADY 
COMMON  ENGBAL 

COMMON  FODTNOK  16)  iFREEWAT 
COMMON  HOLDCAP 

COMMON  IDATE(372),iSN0W,lTABLE 

COMMON  KOUNT 

COMMON   LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSHEOV 

COMMON  PACKTEM.PARTICE.PASTINT,PL0TOBS,PLOTWE,PRECIP,PREWE0V 
COMMON  RADIN,RADLWN,RADSWN,REFLECT 

COMMON   SlMTEMl (3) ,SIMTEH2,SIMTEM3,SN0MELT,SO6SE0V( 372) , 
1  SPRECIP(372) ,SPREQV(372) ,SUBTITL(B) 

COMMON  TC0EFF,TEMPMAX,TEMPMIN,THRSHLD,TITL£(8),T0TPREC 
COMMON  USEMEAN 
COMMON  XHAX 
INTEGER  ACTDATE 
INTEGER  DATE. DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER  PASTINT.PLOTOBS.PLOTHE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

C0MM0N/C0NVERT/FIVE9TH.THIRTY2 

C  COMPUTE   THE   AMOUNT   OF   PRECIPITATION  OCCURRING  AS   SNOW  BY 

C   AMOUNT   SNOW  =   P  *    (1.0  -   B/A).  WHERE 

C   p  =   PRECIPITATION   IN  INCHES 

C   B  =  DAILY  MAXIMUM  TEMPERATURE  -  BASE  TEMPERATURE   (DEGREES  F) 

C   A  =  DAILY  MAXIMUM  TEMPERATURE  -  MINIMUM  TEMPERATURE   (DEGREES  F) 

B  =  TEMPMAX  -  BASTEMF 

A  =   TEMPMAX  -  TEMPMIN 

AMTSNOW  =  PRECIP  *   (1.0  -  (B/A)) 

C  NOW  COMPUTE  THE   AVERAGE   TEMPERATURES   (DEGREES   C)   WHICH  PRODUCE 

C  SNOW  AND  RAIN 

TFORSNO  =   (((TEMPMIN  +   BASTEMF)    *  0.5)   -  THIRTY2)   *  FIVE9TH 

TFORAIN  =    (((TEMPMAX  +   BASTEMF)    *  0.5)   -  THIRTY2)   *  FIVE9TH 
C  COMPUTE   THE   EFFECT  OF  THE   SNOW  ON  THE  SNOWPACK 

CALL   SNOWED   ( TFORSNO, AMTSNOW ) 

C  COMPUTE   THE   EFFECT  OF   THAT   PORTION  OF  THE   PRECIPITATION  OCCURRING 

C   AS  RAIN  ON  THE  SNOWPACK 

CALL   RAINED   ( TFORA I N . PRECI P-AMTSNOW ) 

RETURN 

END 


Subroutine  RADBAL 


Subroutine  PLOTTER 


SUBROUTINE  PLOTTER 

C  THIS   SUBROUTINE  CONTROLS  THE   PLOTTING  OF  THE  MODEL  OUTPUT 

COMMON   ACTDATE. ACTUAL(2l) . AVETEMC , AVETEMF 
COMMON  BASTEMF 

COMMON  CALAIR,CALDEF,CALORIE,CALSNOW,COVDEN 
COMMON  DATE.DATESO)  .DEN, DENSITY, DREADY 
COMMON  ENGBAL 

COMMON   FOOTNOT ( 16) ,FREEWAT 
COMMON  HOLDCAP 

COMMON  IDATE(372),ISN0W,ITABLE 

COMMON  KOUNT 

COMMON  LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON   PACKTEM, PART  ICE .PAST  INT. PLOTOBS.PLOTWE, PRECI P. PREWEQV 
COMMON  RAO  I N.RADLWN.RADSWN. REFLECT 

COMMON   SlMTEMl (3) , S I MTEM2 , S I MTEM3 . SNOMELT . SOBSEQV ( 372 ) . 
1  SPRECIP(372) ,SPRE0V(372) .SUBTITL(8) 

COMMON  TCOEFF. TEMPMAX. TEMPMIN,THRSHLO.TITLE( 8). TOTPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE. DATES. DREADY 
INTEGER  FOOTNOT 

INTEGER  PASTINT,PLOTOBS,PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 
DIMENSION  LEGENDIB) 

DATA  LEGEND/  80H1   =   PREDICTED  WATER  EQUIVALENT,    2  =  OBSERVED  WATER 
1   EOUIVALENT,   3  =  PRECIPITATION/ 

C  TURN  OFF  THE   AUTOMATIC   PAGE  EJECT 

WRITE  (6,910) 
910  FORMAT(«0*) 

C  INITIALIZE   THE   PLOT  ROUTINE 

CALL   TSPLOT    ( 3 , 0 , 0 ,0 ,0 ,0 ,0 , XMAX ,0. U , I  DATE ( 1 ) , 7,-1 ) 

C  PLOT  ONE   LINE   FOR   EACH  INTERVAL 

DO   10   I    =    1, KOUNT 
10  CALL   TSPLOT    ( 3 , SPRE QV ( I ) . SOBSEQV ( I ) . S PRECI P ( I ) . 0 . 0. 0. 0 . 0. 0. XMAX , 
1   0.0, IDATE ( I ) ,7, 1 ) 

C  WRITE   THE   TITLE   AND  OTHER   IDENTIFYING  INFORMATION 

WRITE    (6,920)    T I TLE , SUBT I TL . TCOEFF . COVDEN. LEGEND. FOOTNOT 
920  F0RHAT(1H026XBA10/27X8A10/*  TRANSM I  S  I V I  TY  COEFFICIENT  =*F',.2,T,X 

1  *COVER  DENSITY  ^'F*. 2/ IHO , 1 9X»PL0T   LEGEND  —  *B A 10// IHO 1 3A 10/ IX 

2  3A10) 

C  TURN   THE   AUTOMATIC   PAGE   EJECT  SACK  ON 

WRITE  (6.930) 
930  FORMAT(*R*) 
RETURN 
END 


SUBROUTINE  RADBAL 

C  THIS   SUBROUTINE  COMPUTES  THE   RADIATION   BALANCE  AND  TRANSFERS 

C   CONTROL  TO  THE  DIFFUSION  MODEL    IF    IT   IS  NEEDED 

C  

C  DICTIONARY 

C  

C  SNOCAN  -   THE   LONGWAVE   RADIATION   BALANCE   BETWEEN   THE   SNOW  AND  THE 

C  CANOPY 

C  SNOSKY  -  THE   LONGWAVE  RADIATION   BALANCE   BETWEEN   THE   SNOW   AND  THE 

C  SKY 

COMMON   ACTDATE. ACTUAL (21  I . AVETEMC. AVETEMF 

COMMON  BASTEMF 

COMMON  CALAIR.CALDEF.CALORIE.CALSNOW.COVDEN 
COMMON   DATE .DATES (3) , DEN . DENS  I T Y . DR EADY 
COMMON  ENGBAL 

COMMON  FOOTNOT ( 15 ) ,FREEWAT 
COMMON  HOLDCAP 

COMMON    IDATE(372) , ISNOW.ITABLE 

COMMON  KOUNT 

COMMON  LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON   P ACKTEM, PART  ICE, PAST  I  NT, PLOTOBS, PLOT WE, PRECI P. PREWEQV 
COMMON  RAD  I N.RADLWN.RADSWN. REFLECT 

COMMON   SlMTEMl (3) . S I MTEH2 . S I MTEM3 . SNOMELT . SOBSEQV ( 372 ) . 
1  SPRECIP(372) .SPftEQV(372) ,SUBTITL(B) 

COMMON   TCOEFF. TEMPMAX, TEMPMIN.THRSHLD, TITLE (8). TOTPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER   DATE, DATES. DREADY 
INTEGER  FOOTNOT 

INTEGER   PASTINT, PLOTOBS.PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

C0MM0N/C0NVERT/FIVE9TH,THIRTY2 

C  COMPUTE   THE  CALORIC   INPUT  FROM  NET   SHORT  WAVE  RADIATION  AS  A 

C   FUNCTION  OF   THE  REFLECTIVITY 

RADSWN  =  RAOIN  *    (1.0  -  REFLECT)   *  TCOEFF 

C  IF   THE   PRECIP  WAS   SNOW,    THE   NET  LONG  WAVE   RADIATION  BALANCE  IS 

C  ASSUMED  TO  BE  ZERO 

IF(ISNOW)  20,20,10 
10  RADLWN   =  0.0 

GO  TO  50 

C  TO  COMPUTE   THE   LONG  WAVE  RADIATION  COMPONENTS,   CONVERT  THE  AIR 

C   AND  SNOW   TEMPERATURES   TO   POTENTIAL  CALORIES  BY  THE   STEFAN  - 

C   BLOTZMANN  FUNCTION,   CALORIES  =   S   »   (T   **  4),  WHERE 

C   s  =   1.17E-7  CAL/( (CM*»2) ( DEGREES  KELVIN)**^),  AND 

C   T  =  ABSOLUTE   TEMPERATURE    (DEGREES  KELVIN) 

20  CALAIR   =   1.17E-7   *    ((AVETEMC  +   273.161   **  ) 
USE  =  AVETEMC 

C  IF   THE   SNOWPACK   IS   ISOTHERMAL,   USE  THE   MINIMUM  TEMPERATURE  FDR 

C   COMPUTING   THE   BACK  RADIATION 

IF(CALOEF.EQ.O.O)   USE   =    (TEMPMIN  -  THIRTY2)    *  FIVE9TH 

C  UNDER  NO  CIRCUMSTANCES   HAY  THE   TEMPERATURE   FOR  COMPUTING  THE  BACK 

C   RADIATION   BE   GREATER  THAN  ZERO 

IFIUSE.GT.O.OI   USE   =  0.0 

CALSNOW  =   1.17E-7  «    ((USE  *   273.15)    *»  <») 

C  COMPUTE   THE   LONG  WAVE   RADIATION   COMPONENTS   AS  A   FUNCTION  OF  THE 

C   FIRST,   DETERMINE  WHETHER  THE   SKIES   ARE  CLEAR  OR  CLOUDY 

IF(PRECIP)  30,30,40 

C  WITH  CLEAR   SKIES,   THE  DOWNWARD  LONGWAVE  RADIATION  COEFFICIENT  IS 

C  .757    (RUNOFF   FROM   SNOWMELT,    EH  1 1 10-2- 1406 ,   US  ARMY  CORPS  OF 

C   ENGINEERS,    1960,    PAGE  7) 

30  SNOSKY   =    (1.0  -  COVOEN)    *    ((0.757   *  CALAIR)   -  CALSNOW) 

C  THE   DOWNWARD  LONGWAVE   RADIATION  COEFFICIENT    IS   l.O   BENEATH  THE 

C   FOREST  CANOPY   (OR   BENEATH  CLOUDY  SKIES) 

SNOCAN  =  COVDEN  *   (CALAIR  -  CALSNOW) 

RADLWN   =   SNOCAN   *  SNOSKY 

GO   TO  50 

C  WITH  CLOUDY  SKIES,   WHEN  THE   DOWNWARD  LONGWAVE  RADIATION  COEFFI- 

C   CIENT   IS   1.0   INSTEAD  OF   .757,   THE   ABOVE  THREE   EQUATIONS  MAY  BE 

C   REDUCED  ALGEBRAICALLY  TO  THE  FOLLOWING  SINGLE  EQUATION 

40  RADLWN  =  CALAIR  -  CALSNOW 

C  COMPUTE   THE   CALORIC    INPUT  OR  LOSS  FROM  THE   NET  EFFECT  OF  SHORT 

C   WAVE   AND  LONG   WAVE  RADIATION 

50  CALORIE   =  RADSWN  +  RADLWN 

C  THE   SNOWPACK  TEMPERATURE  DIFFUSION  MODEL    (LEAF,    1970,    STUDY  PLAN 

C   FS-RM-1502,   NO.   224.      ROCKY   MOUNTAIN  FOREST   AND  RANGE   EXP   STA)  IS 

C   INCORPORATED  TO  CONTROL  THE   SNOWPACK  TEMPERATURE  AND  COLD  CONTENT 

C   DURING  NON-ISOTHERMAL  CONDITIONS.      SEE  NOW    IF  THE   DIFFUSION  MODEL 

C   HAS  BEEN   INITIALIZED.      IF   SO,   GO  THROUGH  THE  LINKING  ROUTINE 

C   WHICH   INTERFACES   THE   DIFFUSION  MODEL   WITH  THE  RADIATION  BALANCE 

IF(DREADY)  60,70,60 
60  CALL   LINK    ( IRETURN) 

IF(IRETURN)  70,70,90 
70   IF(CALORIE)  80.90.100 
80  CALL  CALOSS  (CALORIE) 
90  RETURN 
100  CALL   CALIN  (CALORIE) 
RETURN 
END 


Subroutine  RAINED 


SUBROUTINE   RAINED    ( TFOR Al N . AMTRA I N ) 
-THIS   SUBROUTINE  COMPUTES  THE   EFFECT   OF  RAIN  ON  SNOW 
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AMTRAIN 
CALRAIN 


C  DICTIONARY 

C 
C 
C 
C 
C 
C- 


THE  AMOUNT  OF  PRECIPITATION  OCCURRING  AS  RAIN  (INCHES) 
THE  DEPLETION  OF  THE   TOTAL  CALORIE  DEFICIT  BY  THIS  RAIN 
(CALORIES) 

TFORAIN  -  THE  TEMPERATURE   FOR  COMPUTING  THE   DEPLETION  OF  THE  TOTAL 
CALORIE  DEFICIT  CAUSED  BY  THIS  RAIN   (DEGREES  C) 


COMMON  ACTDATE,ACTUAL(21) , iVETEMC , AVETEMF 
COMMON  8ASTEHF 

COMMON  CALAIR,CALD£F,CALORIE,CALSNOW.CO«DEN 
COMMON  DATE,DATES(3 ) , DEN, DENSITY , DREADY 
COMMON  ENGBAL 

COMMON   F0QTN0T(16) .FREEWAT 
COMMON  HOLDCAP 

COMMON   I0ATE(372) .ISNOW.ITABLE 
COMMON  KOUNT 
COMMON   LASTUSD, LINES 
COMMON  NEXTACT 
COMMON  DBSWEQV 

COMMON  PACKTEM,PARTICE,PASTINT,PLOTOBS,PL0TME.PRECIP.PREHEaV 
COMMON  RAD  I N, R AD LWN.R ADS WN, REFLECT 

COMMON  SIMTEMl 13) , S I MTEM2 , S I MTEM3 , S NOMELT , SOBSEQV ( 372 ) , 
1  SPRECIP(372),SPRE0V(372),SUBTITL(8) 
COMMON  TC0EFF,TEMPMAX,TEMPMIN,THRSHLD,TITLE(8),T0TPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER   DATE. DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER  PASTIST.PLOTOBS.PLOTHE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

C  ADD  THIS  AMOUNT  OF  PRECIPITATION  TO  THE  PREDICTED  WATER  EQUIVALENT 

PREHEQV  =  PREWEOV  -i-  AMTRAIN 

C  SEE   IF   THERE   IS   A  CALORIE  DEFICIT   IN  THE  PACK 

IFfCALDEF.LE.O.O)   GO  TO  50 
-COMPUTE  THE  AMOUNT  OF   RAIN  AT  THIS   TEMPERATURE  THAT   IS  NEEDED  TO 
-  WIPE   OUT  THE   DEFICIT  AND  COMPARE   IT  WITH  THE  ACTUAL  AMOUNT 
AMTNEED  =  CALDEF/ ( ( 80- 0  +  TFORAIN)   «  2.54) 
COMPARE   =  AMTRAI'J  -  AMTNEED 
IF(COMPARE)  20,10,40 
-THERE   HAS   JUST  ENOUGH  TO  WIPE  OUT  THE  DEFICIT 
CALDEF  =  0.0 
PACKTEH  =  0.0 
GO  TO  30 

-THERE   WAS  NOT  ENOUGH  TO  WIPE   IT   OUT   COMPLETELY.     JUST  DEPLETE 
THE  DEFICIT 

CALDEF  =  CALDEF  -    ((80.0  +  TFORAIN)    *  AMTRAIN  *  2.54) 
PACKTEH  =  -CALDEF/ I PREWEQV*1.27) 

ADD  ALL  THE  RAIN  TO  THE   PACK  AS    ICE  AND  GET  THE  NEW  HOLDING 

CAPACITY 
PARTICE  =  PARTICE  +  AMTRAIN 
HOLDCAP  =  0.04  *  PARTICE 
RETURN 

C  THERE   WAS  MORE  THAN  ENOUGH  TO  WIPE  OUT  THE  DEFICIT.     ADD  THE 

C   FROZEN  PART  TO  THE   ICE  AND  GET  THE  NEW   HOLDING  CAPACITY 

CALDEF  =  0.0 
PACKTEM  =  0.0 

PARTICE  =  PARTICE  +  AMTNEED 
HOLDCAP  =  0.04  *  PARTICE 
-THE  AMOUNT  OF   RAIN  NOT   FROZEN   IS   FREE  WATER   AND  CONTRIBUTES 


C  

C  

20 

C  

C  

30 


  CALORIC   INPUT  TO  THE  PACK 

FREEWAT  =  COMPARE 

CALL  CALIN    (TFORAIN  »  COMPARE  »  2.54) 
RETURN 

C  ALL  OF   THE   RAIN   IS   ADDED  TO  THE   FREE  WATER  AND  CONTRIBUTES  CAL3RIC 

C  INPUT  TO  THE  PACK 

50  FREEWAT  =  FREEWAT  +  AMTRAIN 

CALL  CALIN    ITFORAIN  «  AMTRAIN  *  2.54) 

RETURN 

END 


Subroutine  RDPACK 


SUBROUTINE  RDPACK 

C  THIS   SUBROUTINE  READS   THE   ACTUAL  OR  OBSERVED  SNOWPACK  TEMPERATURE 

C   CARD  AND   INITIALIZES  THE   DIFFUSION   MODEL    (SUBROUTINE  DIFMOD).  IT 

C   IS  CALLED   WHENEVER  COL  70  OF  AN   INPUT   CARD   IS  NOT   BLANK  OR  ZERO 

COMMON  ACTDATE, ACTUAL(2I) , AVETEMC , AV ETEMF 

COMMON  BASTEMF 

COMMON  CALAI R, CALDEF .CALORIE ,CALSNOW,COVDEN 
COMMON  DATE,DATES(3) , DEN , DENS ITY , DREADY 
COMMON  ENGBAL 

COMMON  F00TN0T(16), FREEWAT 
COMMON  HOLDCAP 

COMMON   IDATE(372) , ISNOW, ITABLE 

COMMON  KOUNT 

COMMON  LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON   PACKTEM,  PARTICE, PAST  I  NT, PLOTOBS,PLOTWE,PREC I P,PREWEQV 
COMMON  RADIN.RADLWN.RADSWN. REFLECT 

COMMON   SIMTEMl (3) , SI MTEM2 , SI MTEM3 , SNOMELT , SOBSEQV ( 372 ) , 
1  SPRECIP(372) ,SPRE0V(372) ,SUBTITL(8) 

COMMON  TC0EFF,TEMPMAX,TEMPMIN,THRSHL0,TITLE(8),T0TPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE, DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER  PASTINT.PLDTOBS.PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 


INTEGER  USEMEAN 

READ    (5,910)    I, ACTDATE, ACTUAL 
910  F0RMAT(A10,1XI6,21F3.1 ) 

C  BE   SURE  THIS   IS  AN  OBSERVED  PACK  TEMPERATURE  CARD  AND  THAT  THE 

C   DATE   IS  THE   SAME   AS  THE  CARD  JUST   READ  BY   SUBROUTINE  READER.  IF 

C  NOT.    ABORT  THE  JOB 

IFd.EQ.lOHOBS   PACK  T)   GO  TO  20 
10  PRINT  920, DATE, NEXTACT 
920  F0RMAT(*1THE   INPUT   CARD  DATED*I7,»   INDICATED  THAT  AN  OBSERVED  PACK 
1   TEMPERATURE  CARD   WAS  NEXT    (COL  70  CONTAINED   A*I2,*.*  /*OBUT  EITHE 
2R   COLUMNS   1  -  10  OF  THE   NEXT  CARD  DO  NOT  CONTAIN  -OBS  PACK   T-  OR  T 
3HE  DATES  ARE  NOT  THE  SAME*) 
CALL   SYSTEM   (999,12H-J0B  ABORTED) 
20   IF(DATE  -  ACTDATE)  10.30,10 

C  GET  THE   ACTUAL   AVERAGE  SNOWPACK  TEMPERATURE    (THE  TEMPERATURES  READ 

C   IN  ARE  ALL   POSITIVE  -  THE  NEGATIVES   WILL  BE   INSERTED  EY  THIS 

C— —  SUBROUTINE) 
30  DUMHYl   =  1.0 

DUMMY2   =  ACTUAL(20) 
DO  60   I   =  1,19 

C  WATCH  FOR   BLANKS  TO  PREVENT  AVERAGING  ZEROES    (A   BLANK   IS  READ  AS  A 

C   NEGATIVE   ZERa,   WHICH  MAY  BE  DETECTED   EY  THE  FOLLOWING  TWO  TESTS) 

IF(ACTUAL( I ) )  50,40,50 
40   IF(SIGN( 1.0,ACTUAL( I ) ) )  70,70,50 
50  nUMMYl   =  DUMMYl  +  1.0 
60  DUKMY2   =  DUHHY2  +  ACTUAL(I) 

C  COMPUTE  THE  AVERAGE 

70  SIMTEM3  =  -  DUMMY2/DUMMY1 

C  FILL   IN  -SIHTEMl-  WITH  THE   OBSERVED  TEMPERATURES.      START  WITH  THE 

C   AVERAGE   AIR  TEMPERATURE 

SIMTEMKl)    =  AMINI    ( AVETEHC  ,  0  .  0  ) 

C  GROUND  TEMPERATURE 

SIMTEM1(3)    =  -  ACTUAL(20) 

C  FIND   THE   MEDIAN  DEPTH  OF  THE   PACK  FROM  THE  WATER  EQUIVALENT 

WEOV  =   PREWEOV  +  PRECIP 

Q  IF  THE   DENSITY   WAS  NOT  SPECIFIED  ON  THE   INPUT  CARD,   COMPUTE   IT  AS 

C   A  FUNCTION  OF   THE   WATER  EQUIVALENT.      (THE   FUNCTION  WAS  DERIVED 

C   FROM   1969   OBSERVED  CONDITIONS  ON  THE  ERASER   EXPERIMENTAL  FOREST) 

IF(DEN)  90,90,80 
80  DENSITY  =  DEN 
GO  TO  100 

90  DENSITY  =   (EXP((0.0179  *  WEOV]    +  3.02))/100.0 

C  DEPTH  =  WEQV/DENSITY 

C  MEDIAN  =  DEPTH/2 

100  DUMMYl   =  WEQV/(DENSITY  +  DENSITY) 

C  INTERPOLATE   BETWEEN  THE   ACTUAL  TEMPERATURES   ON   EITHER   SIDE  OF  THIS 

C   DEPTH  TO  OBTAIN  THE  TEMPERATURE.      (THE  SUBSCRIPT  Of  THE  TEMP 

C   ABOVE   IT   IS  FOUND   BY  TRUNCATING  THE   QUOTIENT  OF  THIS  DEPTH 

C   DIVIDED  BY  THE   DISTANCE   BETWEEN  MEASUREMENTS.      THE  INTERPOLATION 

C  IS  LINEAR) 

SUBSCRP   =  DUMMY1/ACTUAL(21) 

I   =  SUBSCRP 

TRUNCAT  =  I 

SIMTEM1(2)   =  -  ACTUAL(I)   -    ((ACTUAL(I)   -  ACTUAL(I+1))   *    (TRUNCAT  - 
1  SUBSCRP)) 

C  INDICATE   THAT  THE   DIFFUSID'J   MODEL  HAS   BEEN   INITIALIZED  AND  IS 

C  READY  FOR  USE 

DREADY  =  1 

RETURN 

END 


Subroutine  READER 


SUBROUTINE   READER  (lEND) 

C  THIS   SUBROUTINE   READS   A  DATA  CARD  AND  COMPUTES  THE  AVERAGE 

C   TEMPERATURES 

COMMON   ACTDATE, ACTUAL(21) , AVETEMC , AVETEMF 

COMMON  BASTEMF 

COMMON  CALAIR,CALDEF,CALORIE,CALSNOW, COVDEN 
COMMON   DATE,DATES( 3) , DEN, DENSI TY , DREADY 
COMMON  ENGBAL 

COMMON   FD0TN0T(16), FREEWAT 
COMMON  HOLDCAP 

COMMON   IDATE(372) , ISNOW, ITABLE 

COMMON  KOUNT 

COMMON   LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON  PACKTEM, PARTICE , PAST  I  NT , PLOT OBS ■ PLOTWE , PREC I P , PREW EQV 
COMMON  RADIN,RADLWN, RADSWN, REFLECT 

COMMON   SIMTEMl ( 3 ) , S I MTEM2 , S I MTEM3 , SNOMELT , SOBSEQV ( 372 ) , 
1  SPRECIP(372) ,SPREQV(372) ,SUBTITL(8) 

COMMON  TCOEFF, TEMPMAX,TEHPMIN,THRSHLD,TITLE ( 8) .TOTPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE, DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER   PASTINT,PLOTOBS, PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

C0MM0N/C0NVERT/FIVE9TH,THIRTY2 

C  READ  A  CARD   AND  CHECK  FOR  THE  END  OF  FILE 

READ  (5,900)  DATES , RAD  I N, TEMDMAX .TEMPHIN, OBSWEQV, PREC IP , DEN, 
1  NEXTACT 

900  F0RMAT(3I2,F4.0,7X,2F4.1,14X3F5.2,l5XIl) 

IF(E0F(5))  20,20,10 
10   lEND  =  1 

RETURN 
20   lEND  =  0 

DATE  =   DATES(  1  )*10000  +  DATES(2)*100  <■  DATES(3) 
IF(DATE)  70,70,30 
C  COMPUTE  THE   MEAN  TEMPERATURE   IN   FARENHEIT,    THEN  CONVERT   IT  TO 
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C   CENTIGRADE 

30  AVETEMF   =    (TEMPMAX  +   TEMPHINI    *  0.5 

AVETEMC  =    (AVETEMF  -  TH1RTY2)   •  FIVE9TH 

IF(USEMEAN)  50,50,'tO 
40   TEMPMAX  =  AVETEMF 

TEMPMIN  =  AVETEMF 

C  SEE   IF   THE   NEXT  CARD   IS  AN   ACTUAL  PACK  TEMPERATURE  CARD 

50   IF(NEXTACT)  70,70,60 
60  CALL  ROPACK 
70  RETURN 

END 


Subroutine  SNOWED 


SUBROUTINE   SNOWED    ( TFORSNO, AMTSNOW I 

C  THIS   SUBROUTINE  COMPUTES  THE   EFFECTS  OF  A   SNOW   EVENT  ON  THE 

C   SNOWPACK 

C  

C  DICTIONARY 

C  

C  AMTSNOW  -   THE   AMOUNT  OF   PRECIPITATION  OCCURRING  AS   SNOW  (INCHES) 

C  CALSNOW  -  THE  CONTRIBUTION  OF  THIS   SNOW  TO  THE  TOTAL  CALORIE 

C  DEFICIT  (CALORIES) 

C  TFORSNO  -  THE   TEMPERATURE   FOR  COMPUTING  THE   CONTRIBUTION  OF  THIS 

C  SNOW  TO  THE  TOTAL  CALORIE  DEFICIT    (DEGREES  C) 

C  


COMMON  ACTDATE,ACTUAL(21) , AVETEMC, AVETEMF 
COMMON  BASTEMF 

COMMON  CALAIR.CALDEF, CALORIE, CALSNOW, COVDEN 
COMMON   DATE,DATES(3) , DEN, DENSI TY, DREADY 
COMMON  ENGBAL 

COMMON  F00TN0T(16) ,FREEWAT 
COMMON  HOLDCAP 

COMMON  IDATE(372),ISNOW,ITABLE 

COMMON  KOUNT 

COMMON  LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON  PACKTEM,PARTICE,PASTINT,PLOTOBS,PLOTWE,PRECIP,PREWEQV 
COMMON  RADIN,RADLWN,RADSWN, REFLECT 

COMMON   SIMTEMl (3) , S I MTEM2 , S I MTEM3 , SNOMELT , SOBSEOV ( 372 ) , 
1  SPRECIP(372 ) ,SPREQV( 372) ,SUBTITL(B) 

COMMON  TC0EFF,TEMPMAX,TEMPMIN,THRSHLD,TITLE(8) ,TOTPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE, DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER  PASTINT,PLOTOBS,PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 
I  SNOW  =  1 


C  ADD  THIS   AMOUNT  OF   PRECIPITATION  TO  THE  PREDICTED  WATER  EQUIVALENT 

C   AND  GET  THE  NEW  HOLDING  CAPACITY 

PREWEQV  =  PREWEQV  ♦  AMTSNOW 
PARTICE  =  PARTICE  +  AMTSNOW 
HOLDCAP   =  0-04  *  PARTICE 

C  THE  SNOW  FALLliJG  WHEN  THE  TEMPERATURE   IS  BETWEEN  35  AND  32  DEGREES 

C   DOES  NOT  ALTER  THE  CALORIC  DEFICIT 

IF(TFORSNO.GE.O.O)  RETURN 

C  COMPUTE   THE  CALORIE   DEFICIT   FOR   THIS  SNOW   BY  THE  EQUATION 

C   CALORIE  DEFICIT  =  S(I)*DELTA  T*P,  WHERE 

C   S(I)   =  SPECIFIC  HEAT  OF  ICE   (.5  CAL/CM/DEGREES  C), 

C   DELTA   T  =  CHANGE   IN  TEMPERATURE  WITH  RESPECT  TO  FREEZING  (0.0 

C   DEGREES  CENTIGRADE),  AND 

C  p  =  PRECIPITATION  IN  CM   (CONVERSION  FACTOR  =   2.5',  CM/IN). 

C   THEREFORE,   CALORIE  DEFICIT   =  0.5  *   (TFORSNO)   *   (AMTSNOW  »  2.54) 

CALL  CALOSS   (TFORSNO  •   AMTSNOW   »  1.27) 

RETURN 

END 


Subroutine  STORE 

SUBROUTINE  STORE 

C  STORE   THE   INFORMATION  NEEDED  FOR  THE  PLOT 

COMMON  ACTDATE, ACTUAL(21) , AVETEMC , AVETEMF 
COMMON  BASTEMF 

COMMON  CALAIR,CALDEF, CALORIE, CALSNOW, COVDEN 
COMMON  DATE,DATES(3) ,DEN,DENSITY, DREADY 
COMMON  ENGBAL 

COMMON  F00TN0T(16) ,FREEWAT 
COMMON  HOLDCAP 

COMMON  IDATE(372),ISN0W,ITABLE 

COMMON  KOUNT 

COMMON  LASTUSD, LINES 

COMMON  NEXTACT 

COMMON  OBSWEQV 

COMMON  PACKTEH, PART  I CE.PASTI NT, PLOTOBS,PLOTWE,PRECIP, PREWEQV 
COMMON  RAD  I N,RADLWN,RADSWN, REFLECT 

COMMON   SIMTEM1(3) , SIHTEM2 , S I MTEM3 , SNOMELT , SOBSEQV ( 372 ) , 
I  SPRECIP(372) , SPREQV(372) ,SUBTITL(8) 

COMMON  TC0EFF,TEMPMAX,TEMPMIN,THRSHLD,TITLE(8),T0TPREC 
COMMON  USEMEAN 
COMMON  XMAX 
INTEGER  ACTDATE 
INTEGER  DATE, DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER  PAST1NT,PL0T0BS,PL0TWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 
C  KEEP  TRACK  OF  THE  LARGEST  VALUE  FOR  SCALING  THE  PLOT 


XHAX  =   AMAXl    (XMAX, OBSWEQV, PREWEQV) 

C  INCREASE   THE  COUNTER 

KOUNT  =   KOUNT   ♦  1 

IDATE(KOUNT)   =  DATE 

SPRECIP (KOUNT)   =  PRECIP 

SPREQV( KOUNT)    =  PREWEQV 

IF(PLOTOBS)  10,10,20 

C  BY  STORING   A  NUMBER  OUTSIDE   THE   LIMITS  OF  THE  PLOT,    IT  WILL  BE 

C  IGNORED 

10  SOBSEOV(KOUNT)   =  -1.0 

RETURN 

C  THE   OBSERVED  WATER  EQUIVALENT   IS   TO  BE  PLOTTED 

20   S06SEQV(K0UNT)    =  OBSWEQV 
RETURN 
END 


Subroutine  TSPLOT 


SUBROUTINE   TSPLOT ( MX ,X 1 , X2 , X3 , X4, X5 , X5, XMAX , XM IN , I T , INT , IN  I T ) 

C  THIS  ROUTINE   DOES  THE   ACTUAL  PLOTTING 

C.   HWANG,   A.C.   HOGGATT       1   JULY,  1962 
C  ARGUMENTS 

C       MX=NUMBER  OF   VARIABLES  TO  BE   PLOTTED,   LESS  THAN  OR   EQUAL   TO  6 

C     Xl=VALUE   ATTACHED  TO  FIRST  VARIABLE.   PLOTTING  SYMBOL  WILL   BE   A  1. 

C     X2=VALUE   ATTCHED  TO  SECOND  VARIABLE.   PLOTTING  SYMBOL   WILL  BE   A  2. 

C  X3 

C  . 

C  . 

C     X6  AND  SO  ON  FOR  XN 

C     XHAX=UPPER   END  OF   ORDINATE  SCALE 

C     XMIN=LOWER   END  OF  ORDINATE  SCALE 

C      IT  =ABCISSA   VALUE.    (I.E.,    T,   FOR  XT) 

C      INT=ABCISSA   LABELLING   INTERVAL.    (I.E.,    EVERY    ICHTH  LINE  OF  PLOT 
C     WILL   BE   LABELLED  WITH  VALUE   OF    IDY  ON  HORIZONTAL  AXIS) 
C      INIT     ^INITIALIZING  PARAMETER,   USED  AS  FOLLOWS. 

C  INIT       =1,   GRAPH  WILL  COMPUTE   AND  PR9NT  ORDINATE,  PLOO 

C     AND  PRINT   FIRST  LINE  OF  GRAPH.    SUBSEQUENT  CALL   WILL   PLOT  AND  PR+NT 
C     A  LINE  OF  GRAPH  ONLY. 

C  INIT     =-1   USED  TO  READY   SUBROUTINE   FOR   PLOTTING  NEW  GRAPH. 

C     SUBROUTINE   DOES  NO  PLOTTING  OR   PRINTING  WITH  THIS   SETTING  OF  INIT. 
C  IF   THE   VALUE  OF   SCALING  PARAMETERS   XHIN  AND/OR   XMAX   DIFFER  FROM 

C  THE   PREVIOUSLY  GIVEN  ONES  WITHOUT   RESETTING  OF   INIT,        TS  PLOT 

C  WILL  RESET   SCALE   ACCORDING   TO  NEW   XMIN  AND/OR   XMAX  AND   PRINT  OUT 

C  NEW  ORDINATE   POINT  VALUES   TO  AGREE  WITH  SCALING  OF  PLOTTED  POINTS 

C  INTERNAL  VARIABLES 

C     CHARS  CONTAINS  BCD  CHARACTERS   USED   AS   PLOTTING  SYMBOL. 
C     PLOT  HOLDS   THE  HOLLERITH   IMAGE   FOR   ONE   LINE  OF  PLOT. 
C     XN.    SCALAR  ARGUMENTS  ARE   STORED   IN   THIS  LINEAR  ARRAY. 

C  DELTA  IS  A  SCALING  PARAMETER,  EQUAL  TO  THE  RANGE  DIVIDED  BY  110,  AND 
C  IS  RECOMPUTED  WHEN  NEW  SCALE   IS   INDICATED   (WHEN   INIT     =1  OR  WHEN 

C  XHAX   OR   XMIN  VALUE  DIFFERING  FROM  PREVIOUS  VALUE    IS  GIVEN). 

C     PASHIN,   PASMAX,   ARE  FOR  REMEMBERING  PREVIOUS  VALUE   OF  XMIN,  XMAX. 

DIMENSION  CHARS(8) ,    PLOT(lll),   0RDPT(6),  XN(6) 

EQUIVALENCE (0RDPT(6) ,PL0T(6) ) 

DATA     CHARS(1)/1H1/,   CHARS ( 2 ) / 1H2/ ,   CHARS ( 3 ) / 1H3/ ,    CHARS ( 4 ) / 1H4/ , 
*CHAR$(5)/1H5/,   CHARS(6)/1H6/,   CHARS(7)/1H  /,   CHARS ( 8 ) / IH. / 
DATA  NCALLS/0/,NDY/0/ 
IF( INIT) 1000, 1000, 9001 
CALL   WAS   TO   INITIALIZE  ONLY.... 
1000  NDY=0 

NCALLS=0 
GO  TO  701 
COMMENCE  BY  PROTECTING  THE   INDEX  MX 
9001  M=MAX0(MX,1) 
M-HIN0(H,6) 
CONJURE   PLOTTING  CHARACTERS 
3  XN(1)=X1 
XN(2)=X2 
XN(3)=X3 
XN(4)=X4 
XN(5)=X5 
XN(6)=X6 

IF(NDY)  15,15,28 
28   IF (PASMAX-XMAX)  45,40,45 
40      IF(PASMIN-XHIN)  45,80,45 
15  NDY=1 

COMPUTE   AND  PRINT  ORDINATE   POINTS,  SCALE. 
45  AN=(XMAX-XMIN)*.2 

ORDPT( 1 )=XMIN 

NCALLS  =  0 

DO   17  1=2,6 
17  ORDPT(I)=ORDPT(I-1)+AN 

WRITE    (6,100)    (ORDPT(N) ,N=1,6) 

DO   705  N=l,2 

WRITE  (6,101) 
705  CONTINUE 

101   FORHAT(  9X1H.,5(21X1H.)) 

WRITE  (6,800) 
800  FORMAT(  9Xlll( IH. ) ) 

100  FORMAT ( 1HI,4(E10.3,12X),E10.3,11X,E10.3) 
DELTA=( XMAX-XMIN)/110. 
CLEAR  PRINTER  LINE. 
80  DO  99  K=2,lll 

PL0T(K)=CHARS(7) 
99  CONTINUE 

PLOT( 1)=CHARS(8) 
CAUSE   X   TO  BE   PUT  ON  THE    INTERVAL  (1,111). 
20  DO  34  N=1,M 

L0C=1.5t (XN(N)-XMIN)/DELTA 
CHECK  FOR   X  WITHIN  THE   INTERVAL    (XMIN, XMAX), 
IF(LOC)  34,34,811 
811    IF(lll-LOC)  34,814,614 
CHARACTERS  FOR   PRINTING   NOW  GET   DROPPED   INTO  PLACE. 
814  PLOT(LOC)=CHARS(N) 


20 


3*  CONTINUE 
CHOOSE   OUTPUT  FORH.AT. 

NCALLS=NCALLS+1 

IFIMOOISCALLS.INT) )  815,2004,815 
CRANK  OUT  ONE   LINE  OF  PLOT. 
ZOO'.  WRITE    (6,700)    IT,  (  PLOT  (  K I  ,  K=l ,  1 1 1 ) 
GO   TO  702 

700  FORMATdH         16 ,  2H.  .  11 1  Al  I 

815  WRITE    (6,6001    ( PLOT ( K I , K=l , 1 1 1 ) 
600  F0RHAT(9H  lllAll 
702  PASMAX=XMAX 
PAS«IN=XHIN 
CODING  COMPLETED. 

701  RETURN 

END 


Subroutine  WRITER 


SUBROUTINE   WRITER  (CODE) 
-THIS  SUBROUTINE  KEEPS  TRACK  OF  THE   PRINTING  DETAILS 
CO«HDN  ACTDATE,ACTUAL{21) , AVETEHC , AVETEMF 
CONKON  BASTEHF 

COMMON  CALAIR,CALDEF,CALaRlE,CALSNOW,COVDEN 
COMMON  DATE,DATES(3I , DEN , DENSITY , DREADY 
COMMON  ENGBAL 

COMMON  FQ0TN0T(16),FREEHAT 
COMMON  HOLDCAP 

COMMON   IDATEI372) ,ISNOW,ITABLE 

COMMON  KOUNT 

COMMON   LASTUSO, LINES 

COMMON  NEXTACT 

COMMON  DBSWEOV 

COMMON  PACKTEM,PARTICE,PASTINT,PLOTOBS,PLOTWE,PRECIP,PREWEOV 
COMMON  RADIN.RADLHN.RADSWN, REFLECT 

CaHHO!   "IIMTEMIO)  ,SIMTEM2,SIMTEM3,SN0MELT,S0BSE0V(372)  , 


1  SPRECIP(372I ,SPRE0V(372) ,SUBTITL(8) 

COMMON  TC0EFF,TEMPMAX,TEMPMIN,THRSHLD,TITLE(8) .TOTPREC 
COMMON  USEMEAN 
COMMON  XnAX 
INTEGER  ACTDATE 
INTEGER   DATE, DATES, DREADY 
INTEGER  FOOTNOT 

INTEGER  PASTINT,PL0T0BS,PLOTWE 
INTEGER  SUBTITL 
INTEGER  TITLE 
INTEGER  USEMEAN 

C  CHECK   THE   LINE  COUNTER 

IF(LINES  -  481  40,10,10 

C  SEE   IF  THIS   IS   STARTING  A  NEW  STATION  -   IF   SO,   BYPASS  THE  FOOTNOTE 

10   IF(LINE:  -   9991  20,30,30 

C  WRITE   THE  FOOTNOTE 

20  WRITE    (6,9101  FOOTNOT 
910  FORMAT( lH0i3A10/lX3A101 

C  HEADINGS 

30  WRITE    (6,9201    T ITLE , SUBTI TL 
920  F0RMAT(1H150X*SN0WMELT   RUNOFF  SIMULATION  K0DEL*/27X8A10/27X8410 ) 

WRITE    (6,930)   TCOEFF , COVDEN 
930   FORMAT(»  TR ANSM I S 1 V ITY  COEFFICIENT  =*F4.2,74X»C0VER   DENSITY  =*F7.2 
1  ) 

WRITE  (6,9401 
940  F0RMAT(*0*32X»TEMPERATURE    (F)        PRECIP  (IM 
lENERGY         SNOHPACK  PREDICTED*/ 
223X*DATE  MAX       WIN       AVE         DAY  ACCUM 

3AL    (CALl      TEMP    (C)      W.E.  (IN1*/1 

LINES  =  0 

40  WRITE    (6,9501    DATES , TEMPMAX , TEMPMIN, AVETENF, PRECIP , TOTPR EC , RADSHN, 
1  RADLWN, ENGBAL , PACK TEM , PRE WEQV 
950  F0RMAT(1H020XI2,2I3,1X3F6.1,2F7.2,3X3F8.1,5XF5.1,5XF6.21 
LINES  =  LINES  ♦  2 
RETURN 
END 

40  -450     634  6000 


NET  RAD  (CALl 


SHORT 


LONG  B 


111 


35 


Sample  Input 


lZ3456789D12i456789Dl234567a901254567890 123*567890123456789012345676901 234567890 


ELEVATION  9800  FT. 

12569  171 
085  PACK  T 

12669  180 

12769  189 

12B69  198 

12969  207 

13069  417 

13169  423 

20169  234 

20269  243 

20369  441 

20469  447 

20569  453 

20669  279 

20769  465 

20869  468 

20969  471 

21069  294 

21169  478 

21269  481 

21369  300 

21469  488 

21569  304 

21669  306 

21769  308 

21869  501 

21969  312 

22069  314 

22169  310 

22269  509 

22369  510 

22469  510 

22569  511 

22659  291 

22769  513 

22869  614 

30169  279 
DBS  PACK  T  ( 

30269  275 

30369  516 

30469  517 

30569  264 

30669  260 

30769  519 

30869  522 

30969  525 

31069  528 

31169  530 

31269  263 

31369  536 

31469  539 

31569  541 

31669  544 

31769  547 

31869  550 

31969  272 

32069  555 

32169  558 

32269  567 

32369  300 

32469  312 

32569  325 

32669  601 

32769  610 

32869  619 

32969  628 

33069  636 

3316?  562 

40169  196 
OBS  PACK  T  < 

40269  562 


ASPECT  SSE,   SLOPE  30  PERCENT 


310  90  -0  6 
I  45  45  45  45  45  45  45  45  45  45 

400  190  -0 

360  170  -0  3 

280     30  -0 

230  -80  -0 

210-180  -0 

200-140  -0 

190-160  -0 

60-110  -0  1 

220-230  -0 

360-170  -O 

350  -60  -0 

300     10  -0 

250     50  -0 

290  -90  -0 

370-120  -0 

370  -80  -O 

410     20  -0 

390     30  -0 

410     50  -0 

360   160  -0 

370  -40  -0 

360  120  -O  2 

310  -40  -0 

400-110  -0 

430     10  -0  2 

250  120  -0  2 

270  -10  -0 

310  -90  -0 

330  -90  -0 

360  -10  -0 

440     30  -0 

470    60  -0  i: 

260  -70  -0  I 

420  -80  -0 

410  -90  -0  1 
31  31   31  31  31   31  31 

380  -90  -0 

370  -90  -0  I 

380     20  -0 

240  -90  -0  II 

320  -80  -0  1< 

200  -90  -0 

200-190  -0 

210-100  -0  I 

330-210  -O  I 

280-130  -0  I 

270-200  -0 

210-140  -O  i 

370-200  -0  I 

420-170  -0  I 

480-150  -O  I 

490     20  -0 

470  120  -0  I 

290  130  -O  3 

4O0  -80  -0  I 

500  -40  -0 

490  140  -0 

290     60  -0 

290     70  -0 

190     50  -0 

360     50  -0 

550  130  -O 

550  250  -0 

550  240  -0 

600  270  -0 

630  320  -0 

340  300  1060  3 
01  01  01  01  01  01 

560  220  -0 


Ji0369  559 
40469  526 
40569  685 
40669  631 
40769  481 
40869  354 
40969  6*9 
41069  645 
41169  481 
41269  156 
41369  481 
41469  515 
41569  384 
41669  ISO 
41769  601 
4IS&9  540 
41969  718 
42069  722 
42169  704 
42269  576 
42369  612 
42469  725 
42569  562 
42669  466 
42769  618 
4Z869  673 
42969  633 
43069  706 
50169  629 
50269  711 
50369  578 
50469  551 
50569  259 
50669  546 
50769  304 
50869  600 
50969  812 
51069  789 
51169  735 
51269  635 
51369  578 
51469  586 
51569  419 
51669  388 
51769  766 
51869  689 
51969  504 
52069  738 
52169  658 
52269  567 
52369  530 
52469  711 
52569  718 
52669  802 
52769  618 
52869  686 
52969  646 
53069  726 
53169  667 
60169  691 
60269  653 
60369  643 
60469  691 
60569  614 
60669  384 
60769  346 
60869  595 
60969  634 
61069  566 
61169  288 
61269  403 
61369  490 
61469  403 
61569  307 
61669  480 
61769  221 
61869  422 
61969  595 


580  300 
440  260 
600  260 
580  330 
380  180 
480  190 
500  180 
510  260 
520  280 
370  290 
370  270 
450  260 
480  170 
440  190 
450  200 
510  210 
540  260 
570  210 
630  310 
600  360 
620  250 
620  300 
330  150 
320  130 
360  140 
540  170 
520  290 
530  300 
610  300 
680  330 
630  320 
460  320 
450  320 
380  320 
360  320 
500  180 
600  220 
590  300 
630  340 
580  340 
620  360 
610  320 
460  320 
390  310 
620  300 
610  380 
700  390 
700  420 
680  360 
650  340  . 
640  370 
630  350 
760  400 
790  440 
810  480 
760  420 
750  4  00 
740  420 
410  300 
500  240 
620  260 
570  340 
680  390 
730  460 
640  360 
500  360 
660  360 
670  360 
670  340 
400  320 
460  320 
470  310 
470  320 
450  340 
500  340 
460  370 
490  370 
600  350 


21 


620&9  634 
62169  662 
62269  480 
62369  365 
i2*69  394 
42569  374 
62669  451 
62769  806 
62869  ei6 
62969  826 
63069  758 


640  340 
650  310 
530  370 
520  420 
520  350 
390  330 
380  320 
640  310 
720  330 
690  320 
690  310 


12345678901234567890123456769012345678901234567890123456769012345678901234567890 

111  25       55  -480     843  3200 

FRASER  EXPERIHEN7AL   FOREST.  COLORADO 
OEAOHORSE  CREEK,   UPPER  NORTH  SLOPE 


ELEVATION   10500  FT. 


ASPECT  NE,    SLOPE   35  PERCENT 


12569 

59 

200  110 

-0 

77 

1  009 

OBS  PACK  T 

012569 

48  48  48 

48 

48 

48 

48 

48  48 

48 

10  60 

12669 

63 

290  210 

-0 

6 

009 

12769 

67 

250  190 

-0 

49 

009 

12869 

71 

170  50 

-0 

5 

009 

12969 

75 

120  -60 

-0 

4 

009 

13069 

149 

100-160 

-0 

0 

009 

13169 

153 

90-120 

-0 

0 

009 

20169 

88 

80-140 

-0 

11 

009 

20269 

92 

-50  -90 

-0 

16 

009 

20369 

165 

110-210 

-0 

0 

009 

20469 

169 

250-150 

-0 

0 

009 

20569 

173 

240  -40 

-0 

0 

009 

20669 

108 

190  30 

-0 

8 

009 

20769 

181 

140  70 

-0 

0 

009 

20869 

185 

180  -70 

-0 

0 

009 

20969 

189 

260-100 

-0 

0 

009 

21069 

120 

260  -60 

-0 

3 

009 

21169 

198 

300  40 

-0 

0 

009 

21269 

202 

280  50 

-0 

0 

009 

21369 

128 

300  70 

-0 

8 

00^ 

21469 

210 

250  180 

-0 

0 

009 

21569 

133 

260  -20 

-0 

6 

009 

21669 

136 

250  140 

-0 

27 

009 

21769 

138 

200  -20 

-0 

4 

009 

21869 

227 

290  -90 

-0 

0 

009 

21969 

143 

320  30 

-0 

34 

009 

22069 

146 

140  140 

-0 

33 

009 

22169 

146 

160  10 

-0 

11 

00'9 

22269 

242 

200  -70 

-0 

0 

009 

22369 

245 

220  -70 

-0 

0 

009 

22469 

249 

250  10 

-0 

0 

009 

22569 

252 

330  50 

-0 

0 

009 

22669 

144 

360  80 

-0 

16 

009 

22769 

258 

150  -50 

-0 

0 

009 

22869 

262 

310  -60 

-0 

0 

009 

30169 

142 

300  -70 

-0 

23 

1  009 

08S  PACK  T 

030169 

35  35  35 

35 

35 

35 

35 

10  60 

30269 

142 

270  -70 

-0 

3 

009 

30369 

271 

260  -70 

-0 

0 

009 

30469 

274 

270  40 

-0 

0 

009 

30569 

141 

130  -70 

-0 

13 

009 

30669 

140 

210  -60 

-0 

20 

009 

30769 

284 

90  -70 

-0 

0 

009 

30869 

289 

90-170 

-0 

0 

009 

30969 

294 

100  -80 

-0 

0 

009 

31069 

299 

220-190 

-0 

0 

009 

31169 

303 

170-110 

-0 

0 

009 

31269 

152 

160-180 

-0 

4 

009 

31369 

313 

100-120 

-0 

0 

009 

31469 

318 

260-180 

-0 

0 

009 

31569 

323 

310-150 

-0 

0 

009 

31669 

328 

370-130 

-0 

0 

009 

31769 

332 

380  40 

-0 

0 

009 

31869 

337 

360  140 

-0 

0 

009 

31969 

169 

180  150 

-0 

39 

009 

32069 

347 

290  -60 

-0 

0 

009 

32169 

352 

390  -20 

-0 

0 

009 

32269 

361 

380  160 

-0 

0 

009 

32369 

194 

180  80 

-0 

4 

009 

32469 

205 

180  90 

-0 

3 

009 

32569 

215 

80  70 

-0 

8 

009 

32669 

399 

250  70 

-0 

0 

009 

32769 

408 

250  80 

-0 

0 

009 

32869 

418 

320  180 

-0 

0 

009 

32969 

427 

330  170 

-0 

0 

009 

33069 

437 

370  220 

-0 

0 

009 

33169 

429 

430  270 

-0 

0 

009 

40169 

150 

310  200 
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